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'sj- ■ Abstract 
^ ■ 

^ I A valence QCD theory is developed to study the valence quark properties of 

Q"^ ■ hadrons. To keep only the valence degrees of freedom, the pair creation through 

the Z graphs is deleted in the connected insertions; whereas, the sea quarks are 
Q^l eliminated in the disconnected insertions. This is achieved with a new "valence 

Q^' QCD" lagrangian where the action in the time direction is modified so that 

O . the particle and antiparticle decouple. It is shown in this valence version of 

QCD that the ratios of isovector to isoscalar matrix elements (e.g. Fa/ Da and 
Fs/Ds ratios) in the nucleon reproduce the SU{6) quark model predictions 
^ . in a lattice QCD calculation. We also consider how the hadron masses are 

^ I affected on the lattice and discover new insights into the origin of dynamical 

mass generation. It is found that, within statistical errors, the nucleon and 
the A become degenerate for the quark masses we have studied (ranging from 
one to four times the strange mass). The vr and p become nearly degenerate 
in this range. It is shown that valence QCD has the C, P, T symmetries. The 
lattice version is reflection positive. It also has the vector and axial symme- 
tries, the latter leads to a modified partially conserved axial Ward identity. As 
a result, the theory has a U{2Nf) symmetry in the particle-antiparticle space. 
Through lattice simulation, it appears that this is dynamically broken down to 
Uq{Np) X Uq{Np). Furthermore, the lattice simulation reveals spin degeneracy 
in the hadron masses and various matrix elements. This leads to an approxi- 
mate Uq{2Np) X Uq{2Np) symmetry which is the basis for the valence quark 
model. In addition, we find that the masses of N, A,/9, 7r,ai, and ao all drop 
precipitously compared to their counterparts in the quenched QCD calculation. 
This is interpreted as due to the disapperance of the 'constituent' quark mass 



which is dynamically generated through tadpole diagrams. The origin of the 
hyper-fine splitting in the baryon is largely attibuted to the Goldstone boson 
exchanges between the quarks. Both of these are the consequences of lacking 
chiral symmetry in valence QCD. We discuss its implication on the models of 
hadrons. 



PACS numbers: 12.38.Gc, ll.15.Ha, 12.40.Aa, 11.30. Rd 
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1 Introduction 



In addition to its classification scheme, the quark model is, by and large, quite suc- 
cessful in delineating the spectrum, structure, and decays of mesons and baryons. 
One often wonders what the nature of the approximation is, especially in view of the 
advent of quantum chromodynamics (QCD) which is believed to be, after all, the 
fundamental theory of the quarks and gluons. In order to address this question, we 
need to understand first where the quark model is successful and where it fails. 

To begin with, we need to define what we mean by the quark model. We consider 
the simplest approach which includes the following ingredients: 

• The Fock space is restricted to the valence quarks only, i.e. 3 quarks for the 
baryon and a quark-antiquark pair for the meson. Although there are variations 
by including quark self-energy and so on which go beyond the instantaneous 
interaction and will invoke higher Fock space (e.g. q^q for the baryon and q^(f 
for the meson), we will not consider them here. 

• These valence quarks, be they the dressed constituent quarks or the bare quarks, 
are confined in a potential or a bag. To this zeroth order, the hadron wave func- 
tions involving m, (i, and s quarks and anti-quarks are classified by the flavor- 
spin and spatial coordinates according to the SUq{Q) x SUq{Q) x 0(3) group (for 
brevity, we shall refer it as SU{Q)). The wave functions are totally antisjnn- 
metric in the color space for the baryons and symmetric in the color-anticolor 
combinations for the meson. For example, the S- and D-wave baryons are de- 
scribed by the 56-plets and the P-wave baryons by the 70-plets. Similarly, the 
S- and P-wave mesons are described by the 36-plets 

• The SU (6) symmetry is broken down to SU (3) x SU (2) by the residual inter- 
action between the quarks which is weak compared to the confining potential. 
The degeneracies within the multiplets are lifted by these residual interactions. 

Of course, additional breakings of flavor SU (3) due to the quark masses are 
responsible for the detailed splitting within the octet/decuplet baryon multiplets 
and the meson nonets. 

There are many different versions of the quark model which share these attributes. 
They have been called naive quark model, non-relativistic quark model, constituent 
quark model, bag model, etc. in the literature. Here we shall refer them generically 
as the Valence Quark Model with the defining features of the lowest Fock space 
(or valence Fock space) and the SU (6) flavor-spin symmetry, albeit approximate, as 
their common denominator. In this work, we shall concentrate our discussion on 
the light quark systems where the valence quark picture is less well understood. For 
mesons with heavy quarks, such as the charmoniums and upsilons, the valence picture 
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based on the non-relativistic potential model which fits experiments reasonably well 
is confirmed by the non-relativisitc lattice QCD calculations [0, We shall not 

address them in this study. 

Given the definition of the valence quark model, it is easier to understand where it 
succeeds as well as where it fails. For example, with the one-gluon like exchange po- 
tential as the residual interaction between the non-relativistically confined quarks 
to describe the hyper-fine and fine splittings of the hadron masses, the valence quark 
model is very successful in fitting meson and baryon masses 0, |, ^ |lOl and baryon 
magnetic moments |T1|, |T^. It is also successful in delineating the pattern of 
electro-magnetic semileptonic and nonleptonic weak decays |jl5[, and the 



Okubo-Zweig-Iizuka (OZI) rule ||T6|, etc. Similarly it is true for the MIT bag model 
where the relativistic quarks and anti-quarks are confined in a bag with one-gluon 
exchange interaction fl^, |18|. 

It is worthwhile noting that all these are based on the valence quark picture aided 
with SU{6) X 0(3) group for its flavor-spin and space group. On the other hand, 
there are notable failures. For example, it fails to account for the U{1) anomaly (the 
rj' mass), the proton spin crisis and the irNa term. All these problems are associated 
with large contributions from disconnected insertions involving sea-quarks |[19|, ^ pT 



These are places where the OZI rule fails badly. Consequently, it is natural not to 
expect the valence quark model to work in these cases. 

There are also other places where the valence quark model does not work well. 
They include hadron scatterings, couplings, and form factors which are well described 
by models utilizing the chiral symmetry inherent in QCD. Examples of successful 
approaches based on chiral symmetry include vrvr scattering 22, [2^ 

u 



^ ^ vector domi- 
nance KSRF relation p5[, low-energy irN scatterings |26|, ttN scattering 
up to about 1 GeV with the skyrmion ||28|, nucleon static properties [^, electro- 
magnetic 



7l form factors, ttNN form factor p9|, NN interaction, and Goldberger 



Treiman relation |30|. All these have been worked out quite successfully by the parallel 
developments which explore the chiral symmetry of QCD. These include the a model, 
current algebra, PCAC, chiral perturbation theory, and the more recent developments 
incorporating large Nc QCD pi], such as the skyrmion [E^, P3 and the contracted 



current algebra [34 



The common theme of these models is chiral symmetry which involves meson 
cloud in the baryon and thus, the higher Fock space beyond the valence. This cloud 
degree of freedom is essential in the case of vector dominance of EM form factors, 
the pion cloud for the Goldberger- Treiman relation, and the non-vanishing neutron 
electric form factor. Therefore, it is a challenge to understand why the valence quark 
model 'works' without spontaneously broken chiral symmetry, and where the hyper- 
fine splitting in hadron spectroscopy and the constitute quark mass come from. 

From the above discussion, it is clear that the Fock space beyond the valence is 
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important and we mentioned two degrees of freedom, namely cloud and sea. How 
to relate these degrees of freedom back to QCD unambiguously, how to find out 
their roles in physical quantities, and, more importantly, how to relate them to chiral 
symmetry are the main subjects of this paper. It turns out that chiral symmetry plays 
essential roles in light hadron spectroscopy as well as hadron structure. We find that 
both the 'constitutent' quark mass and the hyper-fine splitting in light baryons are 
more of a consequence of spontaneous chiral symmetry breaking than that of gluons 
and sea quarks. 

In Sec. 2, we will define these dynamical degrees of freedom in the Euclidean 
path integral formalism for the hadronic tensor in deep inelastic scattering. In Sec. 
3, we will demonstrate the effects of these degrees of freedom in hadron structure. 
In Sec. 4, we introduce a valence QCD theory which modifies QCD to suppress 
quark- antiquark pair production. We will also explore the discrete and continuous 
symmetries of valence QCD in Sec. 4. In Sec. 5, we adopt a lattice action for 
valence QCD and prove its reflection posivity and hermiticity. The pion mass, the 
pion decay constant, and the current quark mass from the axial Ward identity are 
used to define the zero quark mass limit on the lattice. In Sec. 6, we calculate various 
ratios of matrix elements to check the SU{6) relations. The nucleon form factors are 
calculated and presented in Sec. 7. We then study hadron spectroscopy in comparison 
with that of QCD to explore the origin of the hyper- fine splitting and the 'costitutent' 
quark mass in Sec. 8. Perhaps the most exciting aspect of valence QCD is a new 
understanding of the origin of dynamical mass generation, something missing in the 
valence quark model and put in by hand via a constituent quark mass. In Sec. 9, we 
compare the symmetry breaking patterns in valence QCD and QCD. Finally, in Sec. 
10, we summarize the lessons learned from the valence QCD and draw an analogy 
between the valence quark model and the nuclear shell model. We will also discuss 
the implication on model building for hadrons. 



2 Quark Dynamical Degrees of Freedom 

We have so far alluded to the meson clouds and sea quarks in additional to the valence 
quarks. They appear in various QCD-inspired hadronic models and effective theories. 
How does one define the valence, the cloud, and the sea quarks unambiguously and in 
a model independent way in QCD? It turns out that the best way of revealing these 
dynamical degrees of freedom is in deep inelastic scattering where the quarks show 
up as the parton densities. 

The deep inelastic scattering of a muon on a nucleon involves the hadronic tensor 
which, being an inclusive reaction, involves all intermediate states 

WM,'') = 7^T.i^^)'^'iPn-p-Q){mM0)\^)HM0)\N)sp^nave.. (1) 
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Since deep inelastic scattering measures the absorptive part of the Compton scatter- 
ing, it is the imaginary part of the forward amphtude and can be expressed as the 
current-current correlation function in the nucleon, i.e. 



It has been shown that the hadronic tensor W^jj{q^, v) can be obtained from 



the Euclidean path-integral formalism where the various quark dynamical degrees 
of freedom are readily and explicitly revealed. In this case, one considers the ratio 
of the four-point function (OAr(t) Ji/(af, t2)<//i(0, ti)OAr(O)) and the two-point function 
(OAr(t — {t2 — ti))OAr(0)), where ONif) is an interpolation field for the nucleon with 
momentum p at Euclidean time t. For example, Ois[{t) can be taken to be the 3 quark 
fields with nucleon quantum numbers, 

On= I rf^xe'P^-"V^'=^(")"(x)((^(")''(x))^C75^('^)"(x)), (3) 



for the proton. 

As both t — t2 » 1/AMn and ti » l/AM^, where AMn is the mass gap 
between the nucleon and the next excitation (i.e. the threshold of a nucleon and a 
pion in the p-wave), the intermediate state contributions will be dominated by the 
nucleon with the Euclidean propagator e~^^^(*~'^*2-ti))_ Hence, 



< O(t)/|fe-?%(f,t2)J.(0,ti)O(0) > 



2A/jv 



< O{t-T)O{0) > 



t-t2» 1/AMn 
ti » 1/AMjv 



J2g-Afiv(t-T) 

1 f d^x 



2MnV 



<N\ f -^e-''^-^J,{xM)U^M)\N >, (4) 



where r = ^2 — ^i, f is the transition matrix element (0|OAr|A^), and V is the 3- volume. 
Inserting intermediate states, Wfj_^{q'^,T) becomes 

W,Af,r) = —^Y.(^nf^'(Pn-P + Q){mMO)\n){n\JM\^^ 
ZMn V n 

(5) 

We see from Eq. that the time dependence is in the exponential factor 
^-iEn~EM)T ^ To go back to the delta function 6{En — En + u) in Eq. (|1|), one needs 
to carry out the inverse Laplace transform ||3^, |35| 

T/ rc+ioo , 

W^,{q\ y) = - dre'^^W.^iq^ r). (6) 

1 J c—ioo 
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This is basically doing the anti-Wick rotation back to the Minkowski space. 

In the Euclidean path- integral formulation of Wf^^iQ"^ , t) , contributions to the 
four-point function can be classified according to different topologies of the quark 
paths between the source and the sink of the proton. They represent different ways 
the fields in the currents and J^, contract with those in the nucleon interpolation 
operator Oat. This is so because the quark action and the electromagnetic currents 
are both bilinear in quark fields, i.e. of the form ^M\l/, so that the quark number is 
conserved and as a result the quark line does not branch the way a gluon line does. 
Fig. 1(a) and 1(b) represent connected insertions (C.I.) of the currents. Here the 
quark fields from the interpolators Oat contract with the currents such the the quark 
lines flow continuously from t = to t = t. Fig. 1(c), on the other hand, represents a 
disconnected insertion (D.I.) where the quark fields from and self-contract and 
are hence disconnected from the quark paths between the proton source and sink. 
Here, "disconnected" refers only to the quark lines. Of course, quarks surf in the 
background of the gauge field and all quark paths are ultimately connected through 
the gluon lines. 




Figure 1: Quark skeleton diagrams in Euclidean path integral formalism for evaluating 
Wfj_i, from the four-point function defined in Eq. (^. These represent the lowest 
twist contributions to W^u- (a) and (b) are the connected insertions and (c) is the 
disconnected insertion. 
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The infinitely many possible gluon lines and additional quark loops are implicitly 
there in Fig. 1 but are not explicitly drawn. Fig. 1 represents the contributions of the 
class of "handbag" diagrams where the two currents are hooked on the same quark 
line. These are leading twist contributions in deep inelastic scattering. 

The other contractions involving the two currents hooking onto different quark 
lines are represented in Fig. 2. Given a renormalization scale, these are higher 
twist contributions in the Bjorken limit. Although details of the operator product 
expansion will be given elsewhere [^], it is easy to see why these contributions in Fig. 
2 are suppressed in deep inelastic scattering. Since the hadronic tensor W^y{q^,u) 
in Eq. (Q) is the current-current correlation for the forward matrix element, the 
large momentum that is injected through the current at ti must be carried out by 
the current at t2 to avoid flowing to the nucleon sink at t. In the case of Fig. 1, 
the momentum will simply be conducted through the quark propagators between the 
currents. However, in Fig. 2, no connected quark propagators exist between the 
current. Thus, the momentum has to be exchanged through the gluon propagator, 
leading to the 0(1/(5^) and 0{as) suppression. This is the reason the parton density 
interpretation is valid for DIS which involves diagonal quark propagators, since the 
cross terms from the quark Fock space wave functions are suppressed. From now on, 
we will neglect these "cat's ears" diagrams in Fig. 2. 




Figure 2: Quark skeleton diagrams similar to those in Fig. 1, except that the two 
current insertions are on different quark lines. These correspond to higher twist 
contributions to W^y and are suppressed by l/Q"^. 

In the deep inelastic limit where < 0{1/Q'^) (we are using the Minkowski nota- 
tion here), the leading light-cone singularity of the current product (or commutator) 
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gives rise to a free quark propagator between the currents. In the time-ordered dia- 
grams in Fig.l, Fig. 1(a) /1(b) involves only a quark/ ant iquark propagator between 
the currents; whereas, Fig. 1(c) has both quark and antiquark propagators. Hence, 
there are two distinct classes of diagrams where the antiquarks contribute. One comes 
from the D.I. ; the other comes from the C.I.. It is frequently assumed that connected 
insertions involve only "valence" quarks which are responsible for the baryon number. 
This is not true. To define the quark distribution functions more precisely, we shall 
call the antiquark distribution from the D.I. (which are connected to the "valence" 
quark propagators and other quark loops through gluons) the "sea" quark. We shall 
refer to the antiquark in the backward time going quark propagator between ti and 
^2 in Fig. 1(b) as the "cloud" antiquark. On the other hand, the quark in the time 
forward propagator between ^2 and ti in Fig. 1 (a) includes both the valence and the 
cloud quarks. This is because a quark propagator from t = to t = t{t > 0) involves 
both the time forward and backward zigzag motions so that one cannot tell if the 
quark propagator between t2 and ti is due to the valence or the cloud. All one knows 
is that it is a quark propagator. In other words, one needs to consider cloud quarks in 
addition to the valence to account for the production of cloud quark-antiquark pairs 
in a connected fashion (Fig. 1(a)); whereas, the pair production in a disconnected 
fashion is in Fig. 1(c). 

One important point to raise at this stage is that this separation into valence, 
anti-cloud and sea is gauge invariant (i.e. in the path- integral formalism of Eq. (|^), 
no gauge fixing is required) and topologically distinct as far as the quark skeleton 
diagrams in Fig. 1 are concerned. However, the separation depends on the frame of 
the nucleon. It is expected that the parton model acquires its natural interpretation 
in the large momentum frame of the nucleon, i.e. p > q. Consequently, in the large 
momentum frame, the parton density for the u and d antiquarks come from two 
sources. 

q{x) =q^{x) + q^{x), (7) 

where q^i^) is the cloud antiparton distribution from the C. I. in Fig. 1(b) and qg{x) 
denotes the sea antiparton distribution from the D. I. in Fig. 1(c). The strange and 
charm quarks would only contribute in the D. I. in Fig. 1(c). Similarly, the u and d 
quarks have 2 sources, i.e. 

q{x) = qv+c{x) + qs{x), (8) 

where qv+c{x) denoting the valence and cloud quarks and qs{x) denoting the sea 
quark are from Fig. 1(a) and Fig. 1(c) respectively. Upon defining qdx) = qdx) 
(note the subscript "c" denotes the cloud not charm), the valence parton distribution 
is obtained by 

qv{x) = qv+c{x) -q^x), (9) 

and is responsible for the baryon number, i.e. / uv{x)dx = J[u{x) -'u{x)]dx = 2 and 
/ dv{x)dx = J[d{x) — d{x)] = 1 for the proton. 
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It has been shown [^] that the sea partons in Fig. 1(c) cannot give rise to a large 
Gottfried sum rule violation, i.e. Us{x) = ds{x). Whereas, the origin of u{x) ^ d{x) 
can come primarily from the cloud antipartons in Fig. 1(b). 

After the dynamical degrees of freedom are established in deep inelastic scatter- 
ing, we need to address their relevance to the quark model. The quark model is 
designed to delineate hadron properites in the rest frame or at low energies, such as 
hadron masses, decay constants, form factors, electroweak transitions, etc. Unlike the 
hadronic tensor which entails the calculation of 4-point functions as illustrated in Eq. 
(|), these quantities involve 2-point and 3-point functions. The question is where do 
the dynamical degrees of freedom reside in the 3-point functions which describe the 
matrix elements of hadrons? To track the degrees of freedom, we can consider the 
operator product expansion as an illustration. 

Since the momentum transfer |g| and energy transfer u are large in DIS, the 
product of currents in the forward Compton amplitude T^i,(g^, u) can be expanded as 
a series of local operators. The matrix elements of these local quark bilinear operators 
are then related to the moments of the parton distribution. This is relatively easy to 
carry out in the path integral formalism since it treats the fermions as Grassmann 
numbers instead of operators. While the details are not important here and will be 



given elsewhere all we need to know is that the effect of expanding in terms 
of l/Q"^ pinches the separation of the two currents (i.e. ti and t2 in Fig. 1) into 
one space-time point. Thus the topologically distinct contributions to the four-point 
functions extracted from Fig. 1(a), 1(b), and 1(c) are related to the matrix elements 
obtainable from the three-point functions in Fig. 3(a), 3(b), and 3(c) respectively. 
The latter represents matrix elements of the series of the local operators. Notice 
that for any single matrix element related to the quark bilinear operator ^'r\l/, i.e. 
(A^|r|A^), Fig. 3(c) which inherits the sea degree of freedom from Fig. 1(c) is still 
distinct from Fig. 1(a) and 1(b) and continues to be a D. I.. On the other hand. 
Fig. 1(a) and 1(b) are no longer topologically distinct. In fact, they represent the 
same C. I. for the local operator. Therefore, the valence and cloud degrees of freedom 
from Figs. 1(a) and 1(b) are lumped together in the C. I. of 3-point functions. They 
cannot be separated in a single matrix element in contrast to the case for four-point 
functions. 

What we have shown in this section is that for a fiavor-singlet current ^'F^', the 
matrix element (iV|r|iV') has both the C. I. (Figs. 3(a) and 3(b) represent the same 
C. I. in this case) and the D. I. (Fig. 3(c)). While one can study the sea effect 
directly and separately from the D. I., one cannot separately study the valence and 
the cloud in the Z-graphs since both are included in the C. I. and Fig. 3(a) and 
3(b) are topologically indistinguishable. Similarly one can trace the quark degrees of 
freedom in decay constants and hadron masses which are obtainable from the two- 
point functions. The flavor non-singlet hadrons are obtainable from the C. I. depicted 
in Fig. 4(a). 
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Figure 3: Quark skeleton diagrams in the Euclidean path integral formalism consid- 
ered in the evaluation of matrix elements for the sum of local operators from the 
operator product expansion of J^{x)J^{0). (a), (b) and (c) corresponds to the oper- 
ator product expansion from Fig. 1(a), 1(b) and 1(c) respectively. 




i ^ t 

(b) 



Figure 4: (a) The connected insertion for a meson propagator, (b) The disconnected 
insertion part for the flavor-singlet meson propagator which involves the annihilation 
channel. 
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The flavor-singlet meson, e. g. t]', also involves the D. I. in 4(b). The quark 
propagators in the two-point functions in Fig. 4(a) and 4(b) include the valence and 
cloud (Z-graphs) only; they do not involve the sea contribution. The sea effects come 
only through the fermion determinant in this case. 



3 Sea and Cloud Quarks in Hadron Structure 

Now that we have identified the whereabouts of the quark degrees of freedom, we shall 
first study the sea effects in hadron structure in order to examine why the valence 
quark model succeeds and fails in the places we alluded to in the introduction. Since 
the valence and cloud are tangled together in the C. I., we can only infer the effects 
of the cloud indirectly. 



3.1 Matrix Elements and SU{6) Breaking 
3.1.1 Flavor-singlet g'^ and Fa and Da 



Polarized DIS experiments at CERN [^8|, ^ and SLAC have extracted the flavor- 
singlet axial coupling constant Qa from the gi sum rule, based on the fact that the 
latter is related to the forward matrix element of the axial current via the operator 
product expansion Since the axial current is the canonical spin operator, g'^ is 
thus the quark spin content of the nucleon; i.e. gA = Am -|- Ad + As, where the spin 
content Aq{q = u, d, s) is defined in the forward matrix element of the axial current, 
ipslqi'-ytM^bqlps) = 2M7vs^Ag. 

The surprising result the experiments found is a small gA (0.27(10) |39| and 
0.28(16) 1^), much smaller than the expected value of unity from the non-relativistic 
quark model or 0.75 from the SU{Q) relation (i.e. 3/5 of the isovector coupling 



gA = 1.2574). This has attracted a lot of theoretical attention H^] and the ensuing 



confusion was dubbed the "proton spin crisis" . 

Direct lattice calculations of gA from the forward matrix element of the flavor- 



singlet axial current has been carried out [|20|, The details of these calculations 



have been given elsewhere [^0|, We shall concentrate on the results here. One 
word of caution though is that the present calculations are still in the framework 
of the quenched approximation and the systematic errors due to finite volume and 
finite lattice size have not been properly assessed. Nevertheless, the smallness of gA 
is understood from these calculations [|^, ^ . As explained in the preceding section, 
gA is composed of two components, i.e. gA = 9\{C.I.) + g\{D.I.) where g\{C.I.) is 
obtained from the connected insertion in Fig. 5(a) and g\{D.I.) is obtained from the 
disconnected insertion with the axial current in Fig. 5(b). 
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Figure 5: Quark line diagrams of the three-point function in the Euchdean path 
integral formalism for evaluating from the flavor-singlet axial- vector current, (a) 
is the connected insertion which contains the valence and cloud degrees of freedom 
and (b) is the disconnected insertion which contains the quark sea. 

Lattice calculation |2D[ indicates that each of the u, d, and s flavors contributes 
—0.12 lb 0.01 to the D.I. (Fig. 5(b)). This vacuum polarization from the sea quarks 
is largely responsible for bringing the value of g'j^ from g^{C.I.) = 0.62(9) to 0.25 ± 
0.12, in agreement with the experimental value. This is an example where the sea 
contributes substantially and leads to a large breaking in the SU (6) relation. Thus, it 
is understandable that it should come as a surprise to the valence quark model — the 
latter does not have the sea degree of freedom and has simply ignored it by assuming 
the OZI rule. It is interesting to note that although the sea effect is large, Aug, Ads, 
and As from the D. I. are almost identical in the lattice calculation. As a result, the 
sea effects cancel among different flavors and are not reflected in Fa = (Am — As)/2, 
Da = (Am — 2Ad + As)/2, and the Fa/ Da ratio. This can be seen from Table 
1 where results from the lattice calculations are tabulated and compared with the 
experiments and the quark model predictions. Since Aug = Ads = As, the Fa, 
Da, and the Fa/ Da ratio are identical with or without the sea quarks of the D. 
I.. Partly for this reason, the relativistic quark model prediction based on the input 
of gA = 1.257 and the SU{6), indicated in the last column of Table 1, are quite 
reasonable for Fa, Da, and the Fa/ Da ratio. It was not until experiments revealed 
a large negative As which leads to a very small gA that the true shortcomings of the 
quark model were disclosed. 
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Table 1: Axial coupling constants and quark spin contents of proton in compari- 
son with experiments, the non-relativistic quark model (NRQM), and the relativistic 
quark model (RQM). 





C. 1. 


C. 1. + D. I. 


Experiments 


NRQM 


RQM 


g\=Au+Ad+As 


0.62(9) 


0.25(12) 


0.28(16) H, 0.27(10)11 


1 


0.754 


g\=Au-Ad 


1.20(10) [|] 


1 1.20(10) 


1.2573(28) 


5/3 


1.257 


g\=Au+Ad-2As 


0.62(9) 


0.61(13) 


0.579(25) [H] 


1 


0.754 


Au 


0.91(10) 


0.79(11) 


0.82(5)^, 0.82(6)15^ 


4/3 


1.01 


Ad 


-.29(10) 


- 0.42(11) 


- 0.44(5)11, - 0.44(6) [||] 


-1/3 


-0.251 


As 




- 0.12(1) 


-0.10(5)13,-0.10(4) 








Fa={^u-As)/2 


0.45(6) 


0.45(6) 


0.459(8) PI 


2/3 


0.503 


D A={^u~2Ad+As) /2 


0.75(11) 


0.75(11) 
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Figure 6: The ratio Ra in Eq. (p!0| ) is plotted against the quark mass. The • shows 
the full result with both the C. I. and the D. I. The o indicates the C. I. result only. 
The solid line is valence quark model prediction of 3/5. 
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The role of the sea is clear. How about the role of the cloud then? Since its 
contribution to the C. I. of three-point funtions is entangled with the valence, we can 
not separate it out as is done for the sea. To see its effect indirectly, we consider the 
ratio 



as a function of the quark mass. Our results based on quenched 16'^ x 24 lattices with 
/? = 6 for the Wilson k ranging between 0.154 to 0.105 which correspond to strange 
and twice the charm masses are plotted in Fig. 6 as a function of the quark mass 
ma = ln{4:Hc/ 1^ — 3). The dotted line is the valence quark model prediction of 3/5 
for both the non-relativistic and relativistic cases. For heavy quarks (i.e. k > 0.133 
or ma > 0.4 in Fig. 6), we see that the ratio Ra is 3/5 irrespective of whether the 
D. I. is included (shown as • in Fig. 6) or not (C. I. alone is indicated as o). This 
is to be expected because the cloud/sea quarks which are pair produced via the Z- 
graphs/loops are suppressed for non-relativistic quarks by 0{p/mq) or 0{v/c). As 
for light quarks, the full result (C. I. + D. I.) is much smaller than 3/5 largely due to 
the negatively polarized sea contribution in the D. I. (Table 1 lists the results at the 
chiral limit). Even for the C. I. alone, Ra still dips under 3/5. As we shall see later 
this is caused by the cloud quarks. As we remove these cloud quarks in the Z-graphs 
by an approximation, Ra will become much closer to 3/5, the SU{6) limit. 

As a cross check, we note that g\ from the hyperon semi-leptonic decay analy- 
sis is 0.579(25) which is many sigmas smaller than the relativistic quark model 
prediction of 0.754. Since the D. I. of g\ is independent of the sea flavors u, d, and s, 
its sea contribution cancels out. Thus 

g\ = Au + Ad-2As = {Au + Ad){C.I.) + {Au + Ad-2As){D.I.) ~ ^^(C.J.) (11) 

Indeed, the lattice result shows g^{C.I.) = 0.62(9) which is consistent with g\. There- 
fore, one concludes that the smallness of gA is related to the cloud quark effect to be 
unraveled later. 

3.1.2 nNa Term, ss in Nucleon and Fs and Ds 

Another place where the sea plays an even larger role is the nNa term. Like the 
pion mass in the meson sector, the TrNa term is a measure of the explicit chiral 
symmetry breaking in the baryon sector. There is a long-standing and well-known 
puzzle surrounding the irNa term 



where m = (m^ -I- md)/2. The isospin even irN scattering amplitude at the Cheng- 
Dashen point |^ which equals to aT,N{,(f = 2m^) to 0{ml/m%) is determined to be 




g'}^ _ Au + Ad + As _ (Am + Ad){C.I.) + (Am + Ad + As) (£)./.) 
g\~ Au- Ad ~ Au- Ad 



(10) 



(jTrAT = m{N\uu + dd\N), 



(12) 
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around 60 MeV. On the other hand, if one assumes that (iV|ss|A^) = 0, a reason- 
able assumption from the OZI rule, the o"^'^ obtained from the octet baryon masses 
gives only 32 MeV almost a factor two smaller than a.^Mi.'f' = 2m^) extracted 



from the tiN scattering. It has been suggested that the resolution would involve 



a soft scalar form factor [H8| and a large ss content in the nucleon |47, 49 with 



y = 2{N\ss\N)/{N\uu + dd\N) ~ 0.2 - 0.3. Both aspects of the resolution are 



checked quantitatively in the lattice calculations [pi] , pOf - It is found that the scalar 



form factor is indeed very soft in the D. I. which is consistent with a two vr in- 



termediate state in the dispersion approach of the chiral perturbation theory [^. A 



large ss content in the nucleon is also found ^1, 50| so that the conjecture of large 



ratio y to resolve the irNa term discrepancy is verified. 

A large ss content in the nucleon implies an even larger D. I. for the {p\uu + dd\p) . 
We see from Table 2 that the D. I. of {p\uu + dd\p) is in fact ~ 2 times larger than 
the C. I. counterpart. Thus, 2/3 of the irNa term comes from the sea (Table 2). This 
is the largest sea contribution to any physical quantity that we are aware of. As for 
the quark model, {p\uu + dd\p) is predicted to be < 3. The reason is the following: 
the scalar current expanded in the plane-wave basis. 



/ 



/TTl — * — * — * — * 
d^k— y^[h\k, s)b{k, s) + d\k, s)d{k, s)], (13) 



is the sum of the quark and antiquark numbers weighted by the factor m/E. Since 
the valence quark model does not contain antiquarks and m < E, except for the zero 
modes, the scalar matrix element will be less than the valence quark number. Thus, 
we see that the quark model underpredicts the TiNa term by a factor of 3 or more. 

Unlike the case of the axial current matrix element, different flavor contributes 
differently to the D. I. of the scalar matrix element - s contributes less than u and d. 
As a result, the SU{3) antisymmetric and symmetric parameters, 

Fs = {{p\uu\p) - (A^|ss|A^))/2 = {{p\uu\p){C.I.) + {p\uu\p){D.I .) - (A^|ss|A^))/2 
Ds = {{p\uu\p){C.I.) - 2{p\dd\p){C.I.))/2 + {{N\ss\N) - {p\dd\p){D.I.))/2 (14) 

are strongly affected by the large D. I. part. We see from Table 2 that both Ds and 
Fs compare favorably with the phenomenological values obtained from the SU(3) 



breaking pattern of the octet baryon masses with either linear |5^ or quadratic 
mass relations [Q. This agreement is significantly improved from the valence quark 
model which predicts Fs < 1 and Ds = and also those of the C. I. alone [^, 
The latter yields Fs = 0.91(13) and Ds = -0.28(10) which are only half of the 
phenomenological values |51|, This again underscores the importance of the 

sea-quark contributions. 

Next, we address the effect of the cloud quarks in the C. I.. Similar to the ratio 
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Table 2: Scalar contents, a-j^N, Fs, and Ds in comparison with phenomenology and 

quark model (QM). The 17.7 MeV in the last column is determined with the quark 

mass from the lattice calculation. 

C. I. C. I + D. I. Phenomenology QM 



{p\uu + dd\p) 

{p\uu — dd\p) 

{N\ss\N) 

Fs 

Ds 



3.02(9) 

0.63(9) 

1.53(7) 

0.91(13) 

-0.28(10) 

17.8(5) MeV 



8.43(24) 



1.51(12) 
-0.88(28) 
49.7(2.6) MeV 



1.52 IT 



1.81|3| 



-0.52P, g — -0.57|3| 



45 MeV H 



< 3 

< 1 


< 1 


< 17.7 MeV 



Ra in the axial-vector case, we plot the ratio 



9 s 



1=0 



{p\uu + dd\p) {{p\uu + dd\p)){C.L) + {{p\uu + dd\p)){D.L) 



{p\uu — dd\p) 
as a function of the quark mass in Fig. 7. 
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Figure 7: The ratio Rs in Eq. (|15]) is plotted against the quark mass. The • shows 
the full result with both the C. I. and the D. I. The o indicates the C. I. result only. 
The solid line is non-relativistic quark model prediction of 3. 
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The dotted line is the valence quark model prediction of 3 for both the non- 
relativistic and relativistic cases. Again for heavy quarks (i.e. k > 0.133 or ma > 0.4 
in Fig. 7), we see that the ratio Rs is 3 irrespective whether the D. I. is included 
(shown as • in Fig. 7) or not (C. I. alone is indicated as o). As for light quarks, the 
full result (C. I. + D. I.) is much larger than 3 largely due to the large sea contribution 
in the D. I. (Table 2 lists the results at the chiral limit). Even for the C. I. alone, 
Rs still overshoots 3. As we shall see later, this is caused by the cloud quarks. As 
we remove these cloud quarks in the Z-graphs in an approximation, Rs will become 
much closer to 3, the SU{6) limit. 



3.1.3 Magnetic Moments of the Nucleon 

After having established the importance of the sea and cloud effects in the axial 
and scalar matrix elements, one should naturally ask what happens to vector current 
matrix elements, especially the neutron to proton magnetic moment ratio /i„/ fip. How 
much will the sea and cloud affect the ratio and in what way? After all, the fin/f^p 
ratio was well predicted by the valence picture - a celebrated defining success of the 
SU{6) symmetry. 

I has been known for some time that a nontrivial sea-quark contribution to baryon 
magnetic moments is essential to reproducing the experimental moments |55|, 
It turns out that the individual sea contribution of each flavor is not small [^, ^ . A 
direct lattice calculation of the strangeness magnetic form factor is G%j{0) = — 0.36± 
0.20 1^ [[ Although the central value of our lattice result differs in sign from that of 
the SAMPLE experiment which has Gij{Q^ = O.lGeV^) = +0.23±0.37±0.15±0.19 
from the elastic parity- violating electron scattering |5^, they are consistent within 



errors. The u and d contributions are ~ 80% larger, {0){D.I.) = —0.65 ± 0.30. 
However, their net contribution to the proton and neutron magnetic moment 

fi{D.I.) = {2/3GIj{0){D.I.) - l/3Gi,{0){D.I.) - 1/3^^(0))/.^ = -0.097 ± 0.037/i^ 

(16) 

becomes smaller due the cancellation of the quark charges of u, d, and s. This small 
S'f/(6)-breaking sea quark effect is further nullified by the cloud effect. 

This is illustrated in Fig. 8, where the neutron to proton magnetic moment ratio 
is plotted against the quark mass. We see in Fig. 8 that for heavy quarks in the 
charm region (m^a at 0.55 corresponds to rriq ~ 1 GeV), the ratio from the C. I. 
is close to the SU{6) prediction of -2/3. This is quite reasonable, as this is in the 
non-relativistic regime where one expects SU{6) to work well as in the case of qa 
and Qs in the preceeding sections. When the quark mass comes down to the strange 
region (i.e. m^a = 0.07), the ratio becomes less negative. Extrapolated to the chiral 
limit, the ratio is —0.616 ±0.022 which deviates from the SU (6) prediction by ~ 8 %. 

*Note this is defined via the vector current S7^s without the electric charge 
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Figure 8: The ratio of the neutron to proton magnetic moment ratio /^n/z^p is plotted 
against the quark mass. The o indicates the C. I. result only and the • shows the 
full result with both the C. I. and the D. I. The solid line is valence quark model 
prediction of - 2/3 and the dashed line is the experimental result of - 0.685. 



We understand this deviation as mainly due to the cloud quark effect in the Z-graphs; 
this will be verified later. Now as the sea quark contribution fi{D.I.) in Eq. (|1^) is 
added to the C. I., we find that it tends to counter the cloud effect and bring the 
ratio back closer to the SU{6) prediction. At the chiral limit, the fin/f^p ratio then 
comes down to —0.68 ± 0.04 which is quite consistent with the experimental value 
of 0.685. We should point out that the /i„//ip ratio for the full result (•) is more 
negative at the chiral limit compared with those at other m^a has to do with the fact 
that different extrapolations are used for the C. I. and the D. I.. The C. I. employs 
the linear quark mass extrapolation, as do other observables for the C. I. [^, |21], ^^-j 
whereas, the D. I. uses the ^yfn^ dependence for the chiral extrapolation as in 
the scalar case to reflect the non-analytic behavior in chiral loops |5^, From 
this analysis, we see that although the individual GJ^j{0){D.I.),Gf,j{0){D.I.), and 
G\j{0) are large, their net contribution fi{D.I.) in Eq. (0) is much smaller because 
of the partial cancellation due to the quark charges of u, d, and s. The net sea 
contribution is further cancelled by the cloud effect to bring the fin/fJ'p ratio close 
to the experimental value. Barring any known symmetry principle yet to surface, 
this cancellation is probably accidental and in stark contrast with the ivNa term 
and flavor-singlet g'^ where the cloud and sea effects add up to enhance the SU{6) 
breaking |^0|, ^ as detailed in the preceeding sections. 

So it seems that the success of the SU{6) prediction of the neutron to proton 
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magnetic moment ratio is fortuitous due to these two cancellations and as a result it 
has not revealed, and may even has helped conceal, the importance of the cloud and 
sea quark effects over the years. 



3.2 Form Factors and Meson Dominance 

In all the above ratios we considered, i.e. Ra,Rs, and fin/ fJ'p, the 5^7(6) breaking 
due to the cloud in the Z-graphs is at the level of 10 - 20% which is relatively small 
compared with, say, the sea quark effect in Rs- However, its effect is large in the 
nucleon form factors and has been a subject of wide interest. 

It is well known from the e+e~ — > nir experiment that the process is dominated 
by the p meson pole. Extending this p meson dominance of the photon coupling to the 
space-like part of the pion electric form factor F-^i^q"^ < 0), it accounts for 95% of the 
pion charge radius and is the basis of the 'vector dominance' [^]. This is supported 



by the lattice studies of the pion electric form factor ^ It is found that the 
calculated form factor for different quark masses can be well fitted with a monopole 
form 1/(1 + g^/m^), with rrip being the lattice p mass for the corresponding quark 
mass. Similarly, the dipole form of the nucleon electromagnetic and axial form factors 
is interpreted as the product of two monopoles. For example, the isovector part of 
the nucleon Dirac form factor can be written as [|63[ 




Figure 9: The schematic diagram which depicts the photon coupling to the nucleon 
goes through the p meson in a vector dominance picture. 
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Frif) = lime) - Fnen = Ir^'-"^ 



(17) 



to reflect that the dominating process is the photon couphng to the p meson which 
in turn couples to the nucleon as shown in Fig. 9. 

One monopole in {q^) is the p meson propagator, and the other one is QpNNi.'f') = 



fp 2 to parametrize the pNN vertex (see Fig. 9). By the same token, the isovector 
axial ^rm factor with axial meson dominance takes the form 



where QaiNNiq^) is the aiNN vertex and can be parametrized with a monopole form. 
The isovector pseudoscalar form factor should reflect the pion pole for small and 
has the form 



gU^' 



1 - q^/ml g^NNiO) 
where QnNNiq'^) is the irNN form factor. 

1.2 



(19) 
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Figure 10: The isovector axial form factor g\{q^)i the proton electric form factor 
G^^{q^), the strangeness scalar form factor gs{q^)i and the isovector pseudoscalar 
form factor g%{q^) are plotted as a function of - q^ . 
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Figure 11: The quark line diagram for the C. I. which illustrates the meson domi- 
nance picture with different intermediate meson state corresponding to the respective 
probing current. 



Thus one of the major differences of the various form factors of the nucleon is 
reflected in the mass of the meson which dominates the matrix element in the t- 
channel for the specific current. We plot in the following the isovector axial form 
factor g\{q'^)i the proton electric form factor G'^(g^), the strangeness scalar form 
factor gs{<f ) 01' ^^"^ isovector pseudoscalar form factor g'p{q^) in Fig. 10. 

We see that, since g\{q^) and g%{q^) involve only the C. I. and G^^{q^) is dominated 
by the C. I. |5^, their different behaviors in q^ reflect the p, ai and tt propagators in 
the cloud which serve as the intermediate states in the meson dominance picture as 
depicted in Fig. 11. 
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Given that experimentally rua^ = 1230MeV > rup = 769MeV > = 140MeV 
and if one assumes that gpHN^^f)^ QaiNNi.'f) and g-KNNi.'f) have a similar form in g^, 
the fact that g%{q^) falls off faster than G^^{q'^) which in turn falls off faster than g\{q^) 
is then to be expected. We should mention in passing that both G^i^q"^) and g\{q^) 
shown in Fig. 10 from the lattice calculation p4| , |57| agree with the experiments 
within ~ 6%. This is a clear manifestation of the cloud quark effect through the 
meson cloud. We can attempt to define the meson-nucleon vertex by dividing the 
form factors in Eqs. (p!7[), (p!8D, and ([T9| ) by their respective meson propagators. 
These are plotted in Fig. 12. We see that the resulting gpNNiq"^), gaiNNiq"^) and 
g-KNN{q^) extracted this way are much closer to each other than those in Fig. 10. We 
should mention that the monopole fit of gnNNiq"^) gives gTrNN^O) = 12.2 ± 2.3 which 
checks out the Goldberger-Treiman relation [|65|]. 
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Figure 12: The meson-nucleon-nucleon vertices gpNN^q"^), gaiNN^q"^), g-irNN^q^)^ and 
gaoNNiq"^) deduced from the EM form factors, g\{q^) , g'p{q^) , and the isovector scalar 
form factor g's{q^) are plotted as a function of - q^. 

Also plotted in Fig. 10 is the strangeness scalar form factor gs{q^) which is from 
the D. I. (Fig. 3(c)). It is very soft and has been interpreted as due to the KK 
intermediate state as depicted in Fig. 13(a) [^. The D. I. with u or d quarks are 
even softer [|21[ and are consistent with the dispersion analysis of vrvr intermediate 
state in chiral perturbation theory (x PT) |^8|. This appears to be the source of 
the pion and kaon loops in x PT which are responsible for the non-analytic 
contribution of m^^^ or and m\ in hadron masses (see Fig. 13(b)). 
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Figure 13: (a) The quark line diagram which illustrates the KK intermediate state 
which dominates the form factor gs{(f)- (b) The kaon loop diagram in chiral pertur- 
bation theory. 
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Figure 14: The neutron electric form factor G'^{q^) together with the fit to the ex- 
perimental result (solid line). The o indicates the C. I. contribution and the • shows 
the full result with both the C. I. and the D. I. 
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This nonlinear dependence on the quark masses or has been observed promi- 
nently in hadron masses with dynamical fermions in lattice simulations [^. This 
illustrates the sea quark effect in hadron masses and form factors. The neutron 
charge form factor in the strict SU{6) quark model would be identically zero, since 
the positively charged u quark and the negatively charged d quarks have the same 
spatial wave function. Thus, the small positive G'%{q^) signals the effects of the cloud 
and the sea without the contamination of the valence part like in other quantities. 
We present the lattice calculation of G^(g^) ||5^ in Fig. 14 together with the experi- 
mental result. It is seen that both the cloud from the C. I. and the sea from the D. 
I. are positive and their contributions are similar in size. 



3.2.1 Density and Size of Nucleon 

Now, we can look at the time averaged radial density distribution of the nucleon due 
to different current probes. Define the time averaged density distribution as 

P(r) = / rftrf^ge^(^^-^^"*'*)F(g2) = Nj^^ J d'q5{qo)e^'l' F{q') (20) 

where N is the normalization factor so that / d^rp{r) = 1. We plot the pseudoscalar 
density pp{r), the scalar strangeness density p|(r), the electric charge density Pc{f)i 
and the axial current density p^(r) so obtained from gp{q^) , Qsi.'f) ^ Cl^Ei.'f)^ ^i^d g\{q^) 
in Fig. 15. 

We see that pp{r) has the longest range. This is presumable due to the the 
pion cloud which dominates the pseudoscalar channel and has the longest Compton 
wavelength of all hadrons. The next longest is the scalar strangeness density ps's(r) 
which seems to reflect the KK meson intermediate state in Fig. 13(a) which would 
correspond to the kaon loop in chiral perturbation theory as shown in Fig. 13(b). 
Then comes the electric charge density in the proton pdr) which is well known and 
has been frequently used to extract the size of the nucleon. Finally, the one with 
the smallest size is the axial current density pa('") which reflects the small Compton 
wavelength of the ai meson. 

Now what is the size of the nucleon? As seen from Fig. 15, it is in the eyes of 
the beholder. In other words, it depends on what probe is used to measure it. It 
ranges from 3.56(3) fm for the pseudoscalar density, 1.06(9) fm for the strangeness 
density, 0.797(29) fm for the proton charge density, to 0.627(29) fm for the axial 
current density, indeed a large variation. 

We see that even though the clouds in the C. I. do not break the SU (6) symmetry 
as much as the seas in the D. I. for the scalar and axial currents, they afford a large 
variation in hadron form factors and sizes. Short of these meson clouds, the valence 
quark model simply is not capable of delineating the richness of the various form 
factors. A model like the skyrmion, on the other hand, is capable of detailing the 
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Figure 15: The normalized pseudoscalar density pp{r) (in red), the scalar strangeness 
density p|(r) (in yellow), the electric charge density Pc{r) (in green), and the axial 
current density Pa{i^) (in blue) are plotted as a function of the radial distance from 
the center. 



G-T relation ^7j, the meson dominance of the nucleon form factors ||63, negative 



square charge radius of the neutron [27|, etc. All these are achieved via the ingredient 
of the meson clouds. 



4 Valence QCD 

After having examined the roles of the dynamical quark degree of freedom, we come 
back to the original question of what approximation to QCD the valence quark model 
represents. As illustrated in Fig. 5, the sea is only involved in the D. I. part of the 
three-point function and thus can be isolated as far as its contribution is concerned. 
On the other hand, we have stressed in Sec. ^that the cloud and valence contributions 
are lumped in the C. I. in Fig. 5 and cannot be separated a posteriori. Thus to 
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single out the valence effects and compare to the valence quark model requires an 
approximation to QCD. This can be achieved by forcefully eliminating pair creation 
and annihilation by decoupling the quark from the antiquark. In other words, we 
want to eliminate all the Z-graphs such as a typical one illustrated in Fig. 16. 




Figure 16: A typical Z-graph as a diagram in the time-ordered perturbation. 

We introduce valence QCD (VQCD), a theory which is designed to achieve this 
goal. 

First of all, we shall introduce the particle field u and the antiparticle field v in 
lieu of the Dirac field in the valence QCD lagrangian 

CvQCD = -\f^uF^. - u[^^^D, + ^.D + m]u- vf^^D^ + 7 -3 + m]v. (21) 

Comparing with the QCD lagrangian, the valence version has changed the 74 into 
2^^^ for the particle field u and ^^^^^ for the antiparticle field v. We note that the 
u and V fields do not couple. Now we want to prove that the propagator of the u 
field only propagates forward in time and does not zig-zag in the time direction to 
generate particle-antiparticle pairs. The propagator Su{x, y; A) satisfies the equation 

o'/i -I- 1 -> 
- (^^^4 + 1-D + m)Su{x, y; A) = 5{x - y). (22) 

This can be cast in the integral representation with the static propagator S^{x, y; A4) 
as the bare part of the solution The static propagator S^{x,y; A4) satisfies the 
following equation where there is no propagation in the spatial direction 

- (^^^4 + m)St{x, y; A,) = 5{x - y). (23) 
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It is easy to write down the formal solution for S^x, y; 

-m(x4-y4)U"^'^ 



-6{xi - 2/4)e" 
6{x4 - Vi) 



-P 



-m,(j/4-X4)_ 
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P 



m 



5{x-y) 
6{x-y) 



(24) 



where P 



X 

y 



Pe^ ^vi '^'^^^^ ig path-ordered parallel transport factor in the time 



direction. We see that the usual antiparticle propagation in QCD which involves the 
6'(y4 — X4) is now replaced with (5(x4 — 7/4) in the second term. Thus, S^{x, y; A4) is the 
static particle propagator which moves forward in time only. Now the full propagator 
Su{x, y; A) can be represented in an integral equation in terms of S^{x, y; A4) 



y; A) = y; A^) + / d'zS^ix, z; ^4)7 • DS^{z, y; A). 



(25) 



The kernel 7 ■ D is responsible for hopping in the spatial direction. The full solution 
can be obtained by substituting Su with 5*° iteratively leading to a hopping expansion 
series 

S^{x, y- A) = Sl{x, y- A4) + T dz^ [ dhS',{x, z; ^4)7 ■ DS'^{z, y- A^) 



X4 



Z4, 



+ / dz'^ dzi / dh'd'^zSlix, z'- v44)7 ■ DSl{z', z- ^4)7 ■ DSl{z, y- A^) + 



(26) 



It is clear from this expansion that the time integration variables 

^4 ; '^4 ; • • • ^i"^ se- 
quenced between x^ and y^ due to 9 and 5 functions in Eq. (pi]). A typical term in 
the series is shown graphically in Fig. 17. 




Figure 17: A term in the hopping expansion series in Eq. (^) is graphically presented. 

From this we see that there is no time backward propagation in Su{x, y; A). There- 
fore, there is no pair creation or annihilation in the particle propagator, albeit it still 
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propagates forward and backward in the spatial direction. Similarly, one can show 
that the antiparticle propagator Sy{x,y] A) moves only in the time backward direc- 
tion. 

Although there is no pair creation or annihilation in VQCD, there are still quark 
loops in the spatial direction which could lead to non-trivial dynamical effect via the 
fermion determinant. Since we want to emulate the valence quark model, the sea 
quark degree of freedom needs to be removed also. Therefore, we will not include 
the fermion determinant in the calculation. In other words, VQCD is inherently 
"quenched" . 



4.1 Pauli Spinor Representation 

In the process of replacing the fermion field \1/ in QCD by two Dirac spinors u and v 
in valence QCD (VQCD), we seem to have doubled the degrees of freedom. It turns 
out that half of the degrees of freedom in u and v are not dynamical and thus can 
be integrated over. As a result, VQCD can be represented by Pauli spinor fields. To 
prove this, we first look at the particle part of the fermion action from the VQCD 
lagrangian in Eq. (|2TD and write it in terms of the upper and lower components 

D4 + m ia ■ D \ f Ui \ 
-ia ■ D m \ U2 r 



d x[ui{Di + m )ui 



a-JDf 

m 



+ {u2 + ui^ '-)m{u2 + — ^ -ui)]. (27) 

m m 

After changing the field variables Ui Xi; and U2 + ~^^^^^ ui ^1, the action 
becomes 

Su= f d'x[xi{D^ + m - ^^-^)xi + ^imU (28) 
J m 

Since the part has no dynamics and is quadratic, it can be integrated over, 

leaving the particle action represented by the Pauli spinor xi 

Su= I d^xxiiD, + m - ^^±^^)xi. (29) 
J m 

Similarly, the antiparticle action can also be written in terms of the Pauli spinor 

/D^ -\- (J • B 
d^xx2i-Di + m )x2- (30) 
m 
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We can redefine the Dirac spinor as x = ( ) and rewrite the VQCD lagrangian 

as 

^VQCD = —jF^uF^i, - 




It is clear from Eq. (plf ) that particle field Xi and antiparticle field X2 decouple. This 
also proves that Cvqcd in Eq. ( pT] ) does not double the fermion degrees of freedom. 
After integrating out the non-dynamical d.o.f. it has exactly 4 propagting spinors as 
shown in Eq. (PT]). It is worthwhile remarking that the Pauli form of Cyqco in Eq. 



resembles that of the non-relativistic QCD lagrangian after Foldy-Wouthuysen 
transformation. It has the single time derivative like in the Schroedinger action and 
it contains the ^ term, much like the non-relativistic expansion. However, we 
should stress that the Pauli form of VQCD is not a non-relativistic or other expansion. 
Its form is exact. Furthermore, it does not have the spin-orbit, tensor, and Darwin 
terms as in non-relativistic QCD. 



4.2 Discrete symmetry 

Let's explore the symmetries of VQCD and see if there is any change from QCD. 
First we examine the discrete symmetries: the parity, charge, and time reversal. 

The gluon part of the VQCD is the same as in QCD and there is no need to 
modify the transformation of the gluon field. For parity and time reversal, the u and 
V fields transform the same way as in QCD. Thus in VQCD 

where = {—x,X4),x^ = (x, — X4). It is easy to show that the VQCD action 
SvQCD = / d'^x Cvqcd is invariant under the above parity and time reversal transfor- 
mations. 

As for the charge transformation, we need to take into account the fact that u 
and V are particle and antiparticle fields which should be transformed into each other 
under charge transformation. We find that Svqcd is invariant under the following 
charge transformation: 





^-'f^i ) (34) 
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Thus, with the appropriate definition, VQCD satisfies the C, P, and T invariance. 



4.3 Continuous Symmetry — U{2Np) 

Next, we shall address the continuous symmetries. Since the Dirac structure of the 
time derivative is modified in VQCD, it is no longer Lorentz invariant, although 
it is still translational invariant. This should be acceptable for the purpose of our 
study, i.e. low energy hadron physics near the rest frame. After all, the quark model 
is supposed to be an effective theory of low energy and small momentum transfer, 
unlike the parton model which addresses different kinetic regimes. 

Like in QCD, VQCD has the global vector and axial symmetries. It is invariant 
under the U(l) transformation 

u — y e'"u, V — > e'"'v (35) 
This leads to conserved vector currents 

a^j; = 0, V; = 0, (36) 
where the Noether currents associated with these gauge transformations are 

^; = ^(,2i)^, j; = v^^^y. (37) 

Therefore, the particle and antiparticle are separately conserved. This is in con- 
trast to the conserved current = ^Z7^\E' in QCD where only the difference of the 
particle and the antiparticle numbers or the valence number is conserved, i.e. 

7V„ = I d'x^^,^ = I-0f Y.[bliP)UP) - dl{p}ds{p}] (38) 

where b^S and b/d are the creation and annihilation operators of the particle/antiparticle 
in QCD. 

The axial symmetry of VQCD is realized in the 75 transformation 

u — > e'^^'v, V — > e'^^'^'u. (39) 

The lagrangian Cvqcd with m = is invariant under this transformation which 
transposes the u and v terms in the lagrangian. As a result, one has the conserved 
axial currents 

/ 24 + 1 \ 

= zui (40) 
K = ( Y ) 75«. (41) 
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We should point out that there is no Adler-Bell-Jackiw anomaly |]70[ in VQCD. 
This is so because in VQCD there is no quark loop involving the time direction, hence 
there is no triangle diagram to generate the axial anomaly. With m 7^ 0, the axial 
Ward identities are 

df,A^ = 2mui'j5V, (42) 
d^A^^ = 2mvi'y^u. (43) 

It is useful to consider the particle field u and antiparticle field v like two fiavors 
and define ( = ^ ^ ^ ' ^^^^ VQCD lagrangian in Eq. (^) can be written as 

Uqcd = -\f^uF^, - cf^^^^D, + 7 ■ ^ + m]C. (44) 
At the massless limit, VQCD is invariant under the transformation 

C — y e^°^C, C — ' e*"'"^C, C — ' e'^""^'C, ( — > e'^^'^^C (45) 

where t's are the Pauli spinor in the two-component u, v space. The 4 operators 
I , '^3,15^1, and 75r2 are the generators of the U{2) algebra. So, massless VQCD 
has the U{2) vector and axial symmetries in the particle-antiparticle space. For 
degenerate massless Np fiavors, it has the U{2Np) symmetry. It is in contrast to the 
SU{Nf)l X SU{Nf)r X Uv{l) chiral and t/y (1) symmetry of QCD. In VQCD with Np 
fiavors, the charges = / (Px[ifji(t°' /2)u ± tfj^j^lt"' /2)v] do not form a complete 
SU{Nf) X SU{Nf) algebra because the vector and axial current contain different 
fields. This can also be seen from the states. For massless particle, u satisfies the 
Dirac equation 

f^^^D4 + ^-D]u = 0. (46) 



75f satisfies the same equation 

r74 + l 



D4 + 7-^](75^) = 0. (47) 



Therefore, x± = 75"^) is a solution of the Dirac equation in Eq. (|^), but it 

has different particle content, i.e. it is a mixture of particle and antiparticle. As 
a result, x± does not have a definite handedness, it contains both helicity states. 
From this we conclude that massless VQCD does not have SU{Nf)l x SU{Nf)r 
chiral symmetry as in QCD. Instead, it has the vector-axial U{2Nf) in the fiavor and 
particle-antiparticle space. 
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4.4 Zero Quark Mass Limit 




Even though we have explored the axial symmetry of VQCD in Sec. ( [4 .31) at the 
massless limit, there is a concern that the the zero quark mass limit may be singular. 
This can be seen from the Dirac equation in Eq. (^) for a free quark 

(48) 

This leads to two Laplace equations for the upper and lower components of the particle 
field u 

V' ui = 0, v'm2 = 0. (49) 

There are no time derivatives in these constraint equations and thus no dynamics |[71|| . 
Similarly, ones sees that the ^ field in Eq. (^Sf ) is ill-defined for the m = case. 

To address this problem, we consider the following approach. Let's consider the 
fermion part of VQCD lagrangian with a small mixture of antiparticle part in the 
particle action and vice versa 



_74 + l 74-1 , - ril -r'y^-ln , ^4 + 1 , , - Hi 

Cf = -u[ — - — D^ + e — - — Di + m + 'j- D\u-v[ — - — D^ + e — - — Di + m + 'j- D\v 

Zi Zi Zi Zi 

(50) 

We will then let both m and e go to zero. 

Let's first consider the free quark case. In this case, the fermion lagrangian is 

_( d^ + m a-d \ _( ed^^ + m a ■ d \ 
^F = -u\ ^ ^ \u-v{ ^ ^ \v. (51) 

V a ■ a —eo4 + m I \ a ■ a — c'4 + m I 



This involves two time derivatives. The eigenvalues for u are determined from 

-E + m ia ■ p 
ia ■ p eE + m 



1^:1 )=0. (52) 



which are, for small e. 



E = m+—, (53) 
m 

^ m(l — e) , p^, , 

E = ^ ^-(m + — ). (54) 

e m 

Note in Eq. (|53D, the kinetic energy term p^/m is different from fp' /2m in the non- 
relativistic case. Now if we let e approaches zero faster than m, the second branch in 
Eq. (B3) will be decoupled from the physical spectrum. At this limit, a gap between 
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E = p = and -E = cxo is created. This could pose a problem for perturbation 
treatment around this axial symmetry point. 

However, the situation is modified when the quarks are interacting. In this case, 
the Dirac equation for u is. 



-D4 + m a ■ D 1 ( "1 
a ■ D -eD4 + m M ^2 



(55) 



One of the coupled equation from Eq. ( ^5]) is 



(-D4 (l-e)L)4 )U2 = . (56) 

e e e 

If we let m and e approaches zero at the same rate such that m/e = A ^ ^qcd, the 
left-hand side of Eq. (^6]) leads to a constraint equation 



(D^ + a ■ B)u2 + z^ff ■ ^ni = 0. (57) 

The right-hand side leads to two equations, both with linear time dependence 

DiU2 = 0, (58) 
(£>4 - A)m2 = 0. (59) 

Since A ^ ^qcd, the solution from Eq. ( |59D is decoupled from the physical system 
of the hadrons. 

Therefore, the Dirac equation for the interacting massless quark with the e regu- 
lator leads to the following coupled equations, 

D4U1 + id ■ Du2 = 0, (60) 
D4U2 = 0, (61) 

with Eq. ( pTl) as a constraint. This should admit propagating solutions. Similar 
situation exists for v. Thus, we can approach the interacting massless quark case 
with the help of the infrared e regulator. 



5 Lattice VQCD 

In order to solve VQCD, we need a lattice version like in QCD. To this end, we shall 
devise a lattice VQCD action based on Wilson's action in QCD. We shall use the 
following lattice VQCD action 

^VQCD = Sg + Sp + Sp (62) 
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where Sg is Wilson's gauge action which is preserved here and the quark action Sp 
and the antiquark action Sp are 

Sp = '^[u{x)u{x) — k{u{x + a4)(l + '')^Ul{x)u{x) + Uqu{x){1 — '^^u{x)) 

X 

-kY^{u{x + ai){l + -ii)Ul{x)u{x) + u{x){l - ji)u{x + a^))] (63) 

i 

Sp = ^[v{x)v{x) - k{v{x){1 - 74)f/4(x)u(a; + 04) + uov{x){l + 74)1; (x)) 

X 

-i%^{v{x + aj)(l + 'ji)Uj{x)v{x) + v{x){l - ji)v{x + Oj))] (64) 

i 

where uq is the tadpole contribution of the gauge link f/^ which we take to be 
(TrOy/^ \72\. This has the VQCD in Eq. (|2T|) as the classical continuum limit. 



5.1 Reflection Positivity and Hermiticity 

Similar to the continuum case in Sec. ( [4.2| ), lattice VQCD action in Eq. (^) is 
invariant under the corresponding lattice C, P, and T transformations. 

For Euclidean action, it is imperative that it satisfies Osterwalder-Schrader reflec- 



tion positivity in order to allow the Euclidean correlations to be continued back 
to the Minkowski space. We shall follow the derivation for the Wilson action . To 
prove reflection positivity, one needs to show 

{{QF)F) > (65) 

where F is a function of the fields u, u, v, v, and U on the positive time part of the 
lattice and 6 is the time reflection operator. We shall consider the 'link-reflection' 
case where the time reflection is with respect to the t = (0 — ^ 1) link. In this case, B 
is deflned by the transformation 

0Mx,t = '"x,l-t74, 0Mx,t = 74Mx,l-i, (66) 
QVx,t = Vx,l-a4., QVx,t = 74Vx,l-t, (67) 

eU{x,t^;yty) = U\x,l~t„yA-ty). (68) 

Let's illustrate the proof with the u part of the action alone. The proof can be 
extended similarly to include the v fleld. Denoting the fleld variables in the half- 
space with positive time t > 1 hj u^,u^, U~^, and in the other half-space with t < 
by u~,u~, U~ , the VQCD action can be separated into three parts 

SvQCD = S4u+, n+, f/+] + n", U~] + S,[u+, u~], (69) 

where 

S,[u^,u-] = + 74)%,o], (70) 
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is the action which connects S+ and and involves hnks between t = and t = 1. 
Here are have used the temporal gauge. Now since 

QS+[u+,u+,U+] = Sl[Qu+,Qu+,QU+] = S_[u~ ,u~ ,U~], (71) 

the integral in Eq. (^) is then 

{{eF)F) = Z-^ J dUdu+du+e-^+^''^'''^^^^^F[u+,u+,U+] 

Consider the Taylor expansion of the last exponential in Eq. ( |72D 

l-'^E^^i(74 + l)(0^^,i) + - (73) 

The only terms that survive the Grassmann integration are the first two terms and, 
with a diagonal representation of 74, they give semi-positive definite contributions to 
{{QF)F). Extension to include the v field is straightforward and thus the reflection 
positivity for the VQCD action is proved. 

In constructing meson propagators, the usual practice is to first invert the quark 
matrix to obtain the quark propagator from the source to all lattice points, i.e. 
M~^{x, 0). Then the anti-quark propagator which goes backward in time is obtained 
through the hermiticity relation 

M-i^(0,a:)=75M-i(a:,0)75, (74) 

where f indicates the hermitian conjugation in the color and Dirac indices. In VQCD, 
a similar situation exists. In constructing a qq meson, one needs the quark propaga- 
tor M~^{x,0) and and the anti-quark propagator M~^(0,x). It turns out that the 
hermiticity relation 

M-iV0,x)=75M-i(x,0)75, (75) 

still exists so that one can obtain the anti-quark propagator from the quark to con- 
struct a meson propagator as before. 

5.2 Free Quark Propagator 

It is useful to understand the free quark spectrum and its residue at the pole for the 
lattice VQCD and see how different it is from the Wilson and the continuum ones. 
The inverse of the free quark propagator of VQCD in momentum space is 

Sr\p) = 1-^(1- 74) - k{1 + 74)6-^^'^'^ - E[(l + + (1 - l^)e'P% (76) 
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where a is the lattice spacing. We can compute the propagator in discrete time t = rita 

SF{t,p)= r ^SUp)e'^'' (77) 



For rit > and p in the 3 direction, 



Spit^Ps) = - «:e^" + «:(e^'^ - 1)74 - 2^^:73 sin (pga)], (7J 

JD 



where 



A = 1 — 5k — 2k cos(p3a), 
B = A^ - + Ak^ sin^(p3a), 

E is the energy. For small ^30, i.e. p^a <^ 1, 

_ ma + 2 2 / \ 

Ea = ma+ —riPsa) , (80) 

2ma{ma + 1) 

where 

, , 1 — 6/t , , , 

ma = ln(— ), (81) 

is the free quark mass which is the same as in the Wilson case and ma = ^ — 4 is 
the small mass approximation for ma. 

We plot in Fig. 18, the dispersion relation of Ea vs. psa for a range of of ma (ma 
= 0.05, 0.1, 0.5, and 1.0) for both the Valence (dashed lines) and the Wilson case 
(solid lines). We see that for heavy quarks, i.e. ma = 0.5 and 1, the two curves at 
the top are close to each other and they differ at small mass and low momentum. At 
small psa, the behavior of Eq. ( ^0|) holds for the valence case. At ma = ma = 0, 
there is a sigularity at Ea = p^a = 0. For any finite p^a, Ea = ln3 which resembles 
the infinite gap in the massless free quark situation in the continuum (see Eq. (|53D ). 

Finally, we see at zero momentum, the static propagator is 

S{t > 0,p= 0) = -J—e-"^'^^. (82) 

1 — OK Z 

This is the same as in the Wilson case and the wave function normalization factor 
is also the same. To convert lattice matrix elements of local currents with bi- 

1— 6ft _ 

linear quark fields, e.g. ^'(x)r^'(x) to the continuum ones, besides the finite lattice 
renormalization one needs to multiply the factor '^~2k^° = moc"*''^ to take into account 
the finite mass normalization due to the Wilson quark action with tadpole improve- 



ment 1 72, 44 1 . Here ^0 = g~ where Kc is the critical k at which point the pion mass 
is zero and m^a = 1^(2^ — 3) is the tadpole improved definition of the bare quark 



mass in Eq. (pi 



35 



' ' ' I 5 ' ' ' t ' ' ' 1 I 5 ' ' ' t> ' ' ' 2 I 5 ' ' ' 

pa 



Figure 18: The dispersion relations between Ea and p^a for a free quark are compared 
between the Wilson (solid lines) and the lattice VQCD (dashed lines) version for ma 
= 0.05, 0.1, 0.5, and 1.0. 

5.3 Lattice Details 

Before we present our results on VQCD from the Monte Carlo calculation, we should 
mention the lattice specifics and the details of the calculation. 

We use the same gauge configurations which have been used for the study of 



hadron masses, matrix elements and form factors [44, 65, 20, 21, 57] in the quenched 
approximation. This way we keep the scale of the lattice spacing unchanged. These 
quenched gauge configurations were generated on a 16'^ x 24 lattice at /? = 6.0. The 
gauge field was thermalized for 5000 pseudo-heatbath sweeps from a cold start and 
100 configurations separated by at least 1000 sweeps were used. Periodic boundary 
conditions were imposed on the quark fields in the spatial directions. In the time di- 
rection, fixed boundary conditions were imposed on the quarks to provide larger time 
separations than available with periodic boundary conditions. As long as the time 
separations ti and t — ti in Fig. 5(a) are large enough, the form factors should not 
depend on the nucleon interpolation field. All quark propagators in the quenched ap- 
proximation were chosen to originate from lattice time slice 5; the secondary nucleon 
source was fixed at time slice 20 (except for k = 0.154 where the quark propagators 
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from time slice 3 to 22 were used). In the case of VQCD, all quark propagators 
originate from time slice 2 and terminate at time slice 22 for the three-point function 
calculation. We also averaged over the directions of equivalent lattice momenta in 
each configuration; this reduces error bars. 

We have verified that the time separation is sufficient so that there is a plateau 
for the quark bilinear current insertion at time slices ti after the nucleon ground state 
is achieved. The quenched approximation part is done for the lightest quarks with 
K = 0.154, 0.152, 0.148, and 0.140, and g^a^ up to A{27r/Lf. The nucleon masses Mj^a 
for K = 0.154, 0.152, and 0.148 are 0.731(11), 0.883(9), and 1.153(7) respectively. The 
corresponding pion masses m^a are 0.375(4), 0.487(3), and 0.679(3). Extrapolating 
the nucleon and pion masses to the chiral limit where we determine Kc = 0.15672(4) 
and the nucleon mass at the chiral limit to be 0.536(13). Using the nucleon mass 
to set the scale which we believe to be appropriate for studying nucleon properties 
P^ , |65| , pm , the lattice spacing = 1.75(2) GeV is determined. The three k's then 
correspond to quark masses of about 120, 205, and 370 MeV respectively. 

Since we use the same guage configurations for VQCD, the lattice spacing is the 
same as that in the quenched approximation. This is certainly obvious if we choose 
the string tension or the glueball mass to set the scale. Using the physical nucleon 
mass to set the scale in the quenched approximation opens up the question as to 
what extend the fermion determinant effects are implicitly inlcuded. It is shown 



that the quenched approximation can be viewed as including the leading terms in 
the loop expansion of the fermion determinant which are commensurate with the size 
of loops in the gauge action. This leads to a shift in (3 or the coupling constant. 
However, when the infinite volume and continuum limits are taken [^, the scales 



due to hadron masses and the string tension are consistent. Since we are not at the 
infinite volume and continuum limits, the scale set by the nucleon mass differs from 
that by the string tension by ~ 20%. This is the caveat that we should be aware of. 
Nevertheless, whatever scale we decide to choose, the lattice spacing is the same in 
the following VQCD calculation as that in our quenched approximation. 

The determination of Hc which corresponds to the zero quark mass will be dis- 
cussed in the next section. To determine the finite quark mass, we shall use the 
tadpole improved form of the lattice free-quark mass in Eq. (BTI), i.e. 



1 4k 
m„a = ln( 3) = ln(— - 3), (83) 

ZKUn K 



where we have used uq = 



5.4 Pion Decay Constant, Pion Mass and 

Pion decay constant /^r plays an essential role in understanding low-energy chiral 
dynamics and chiral symmetry breaking in QCD. It sets the scale for chiral per- 
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turbation, and relates the Goldstone boson mass to the quark mass through the 
Gell-Mann-Oakes-Renner's relation 



fn^n = -{'mu + nid) {qq) , (84) 

where (qq) is the quark condensate, which is the order parameter for chiral symmetry 
breaking. In VQCD, it is not clear if there is a corresponding relation, or more 
importantly, if the U{2Nf) symmetry is broken to generate Goldstone bosons. We 
can, however, look for clues from the pion decay matrix element with the axial current. 
In QCD, the pion decay constant is defined by 

{0\A^{x)\7r{p)) = zUp^e'P-^. (85) 

Applying axial identity from Eq. (|43|) to the zero-momentum pion state, we obtain 
in VQCD 

{0\d,Al{x)\7r{0)) = -mlU{m^)e-"'-' 

= 2m(0|w75M|7r(0))e-""^*, (86) 

where m is the quark mass. From this, we find 

— ml— ; — -—-= 2m. 87 



It is clear from Eq. (|87|) that as long as the ratio f^(m^)/ (0|i;i75n|7r(0)) does not 
diverge as fast as l/m"^ when the quark mass approaches zero, the pion mass will go 
to zero at massless quark limit. Furthermore, if the pion decay constant is not 
zero at the massless limit, it would signal spontaneous axial symmetry breaking with 
the pion as the Goldstone boson. f„,mT,^ and the current quark mass m from lattice 
VQCD are calculated in the standard way by fitting the following two-point functions 
to 



{El»^75«(x-.()]P(0.0)),^-;_,^- 

X 

{j:mv^^,u{x,t)]P{0,0)) ^-^^ ^-^e~-^\ (89) 

(X:mt)P(0,0))^^-^^^e--*. (90) 

Here P is the pseudoscalar interpolation field ui'j^v and Zj^ = (7r|P|0) is the wave 
function overlap. We use the local current for the axial current in Eq. ( P5| ) for 
the lattice calculation. There are finite lattice renormalizations associated with the 
operators in these matrix elements. We have not calculated them, but we expect the 
multiplicative renormalization constants Za and Zp for the axial and pseudoscalar 
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operators to be of order 1, as in the quenched approximation. Our results presented 
below are subject to this caveat. 

With Wilson-type fermions, one needs to find out where the Kc is for zero quark 
mass. To determine Kc, we plot the dimensionless pion mass, the current quark mass, 
and the pion decay constant in Fig. 19 as a function of m^a = ln(4Kc/fi; — 3) where Kc 
is to be determined from the extrapolation. We see that the pion mass is very linear 
in the range of the quark mass that we calculated. At the same time, the pion decay 
constant behaves like l/m,r in this range. Since from Eq. (|88|) and Eq. (|89[) we 
have the relation 

Uim^)ml = 2mZ^, (91) 
if the linear behavior between m^r and m persists all the way down to zero quark 

Valence 

1-2 I \ 1 1 1 1 1 

current quark ma kOm 
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Figure 19: The dimensionless pion mass rriT^a, the current quark mass ma, and the 
pion decay constant /^a are shown as a function of m^a = ln(4Kc/K — 3). The solid 
lines represent the linear extrapolation with respect to m^a. 

mass and remains constant, then will diverge like 1/m^ or 1/m. Alterna- 
tively, at smaller quark mass than we calculated here, f^^ could conceivably turn 
fiat and in this case m,r will fall off like y/m as in QCD with a constant Z^. Un- 
fortunately, using conjugate gradient to invert the quark matrix, we have encoun- 
tered critical slowing down. The smallest quark mass we run at k = 0.162 already 
takes 5,000 iterations to converge. It is impractical for us to go down any further. 
Short of theoretical guidance and numerical evidence, we extrapolate the pion mass 
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to zero both linearly and quadratically with respect to m^a = ln(4Kc/ n — 3) with 
K = 0.162, 0.1615, 0.1610, 0.1590, 0.1585, 0.1580, 0.1575, 0.1570, 0.1565, 0.1560, 0.155, 
0.154,0.152, and 0.148. We found that = 0.1649(10) (x^ = 0.002 with 14 data 
points) for the linear dependence and Kc = 0.1636(19) (x^ = 0.04 from the three 
largest k's) for the quadratic dependence. We plot in Fig. 20 the quadratic fit of m-„ 
as a funciton of m^a with determined from the linear fit. We see that the Kc 
point from the quadratic fit crosses the abscissa at rriqa = 0.031, however, its error 
bar overlaps with that from the linear fit. Also plotted in Figs. 19 and 20 is the 
current quark mass ma from Eq. (|89| ) as a function of m^a. Extrapolating the quark 
mass linearly with respect to rriqa, we obtain Kc = 0.1642(9) (x^ = 3.5 from the first 
8 k's). We have not used covariance matrix in these extrapolations. We see that the 
current quark mass ma from Eqs. ( [87| ) and (^) crosses the abscissa at m^a = 0.017. 
This is consistent with with that extrapolated from the pion mass, either linearly or 
quadratically. The Kc so obtained overlaps with both of the above two k'^s within 
errors. It is gratifying to know that different definitions of agree. On the other 
hand, it does not differentiate the two scenarios of the pion mass dependence on the 
quark mass. We shall use the linear extrapolaiton with = 0.1649(10) to define zero 
quark mass in this study. Also plotted in Fig. 20 is fT^m^a^ . We see that it is 
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Figure 20: The dimensionless pion mass m^a, f^^m'^a^, and current quark mass ma 
are shown shown as a function of m^a = ln[4Kc/ k — 3) with Kc determined from the 
linearly fit of m^^ with reepect to m^a. 
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quite linear in the range of quark mass that we have considered. This confirms that 
f-jTml oc rrig or equivalently Eq. (pT]), since we have just shown in Fig. 20 that m^a 
and ma are hnearly related. We should stress that we still do not know the behavior 
of the pion mass and pion decay constant when the quark mass is small. But, at 
least we can say that is non-zero (divergent or not) and approaches zero at 
the massless quark limit. This we take to be the evidence that there is a spontaneous 
axial symmetry breaking with the two pions ui'y^v and vi'y^u as the Goldstone bosons 
for each flavor. 



6 SU{6) Relations 

We shall first examine the ratios Ra in Eq. ([10[) , Rs in Eq. (ITsI) and the neutron 
to proton magnetic moment ratio n'^/n^ to check their SU{6) relations as compared 
to those in full QCD, albeit with quenched approximation. In VQCD, the sea quark 
contributions (e.g. Fig. 5(b)) are scrapped. We will only consider the connected 
insertions (e.g. Fig. 4(a) and 5(a)). 



6.1 Ra and Fa/ Da 

The same 100 gauge configurations used for the earlier Ra calculation in Sec. p.l.l| 
are used for the VQCD case. Since in VQCD, there is only the C. I., the Ra ratio 
becomes 

g-X [Au- Ad){C.I.) 

We plot it in Fig. 21 as a function of the dimensionless quark mass niqa (with n 
= 0.162, 0.1615, 0.1610, 0.1590, and 0.1585) together with the corresponding results 
from the quenched calculation presented earlier in Fig. 6 in Sec. p.l.l| . We see that, 
even for light quarks in the strange region (m^a ~ 0.07), it is much closer to the 
valence prediction of 3/5, in contrast to the QCD calculation with C. I. alone. 

This shows that VQCD indeed seems to confirm our expectation of the valence 
quarks behavior, i.e. obeying the SU{6) relation. The deviation from the exact 3/5 
prediction in Fig. 21 reflects the fact that are still spin-spin interaction between the 
valence quarks as evidenced in the a ■ B term in the VQCD action with Pauli spinors 
in Eq. ([29|). Its effects, however, appear to be small. This also confirms our earlier 
assertion that the deviation of the C. I. of Ra in QCD (o in Fig. 21) is largely due 
to the the cloud quark-antiquark pairs. 

With only the C. I., the Fa/ Da ratio is related to VQCD Ra in Eq. (H) 
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Figure 21: The ratio Ra between isoscalar and isovector in VQCD and QCD are 
plotted against the dimensionless quark mass m^a from the strange to the charm 
region. A indicates the VQCD case, o/« indicates the C. I./ Sea + C. I. in the QCD 
case. The dashed hne is the SU{6) prediction of 3/5. 

From Ra = 0.566(11) for the smallest quark mass {n = 0.162), we obtain Fa/Da = 
0.643(4). This is larger than the QCD prediciton of 0.60(2) for the C. I. (see Table 
1) and closer to the valence quark prediction of 2/3. 
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6.2 Rs and Dg/Fs 



Similar results are obtained for the Rs ratio in VQCD. In this case, 



{{p\uu + dd\p)){C.I.) 
{{p\uu-dd\p)){C.L] 



(94) 



Plotted in Fig. 22 are the VQCD results together with those from QCD in Fig. 7 
as a function of the quark mass m^a. Again, we find that the ratios for the light 
quarks are approaching the valence quark prediction of 3. This again confirms that 
the deviation of the C. I. result in QCD is primarily due to the Z-graphs with cloud 
quarks and antiquarks. When they are eliminated in VQCD, Rs becomes close to 
the SU (6) relation. 
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Figure 22: The ratio Rs between isoscalar and isovector scalar charge in VQCD and 
QCD are plotted against the dimensionless quark mass rriqa from the strange to the 
charm region. A indicates the VQCD case and o/« indicates the C. I./ Sea + C. I. 
in the QCD case. The dashed line is the valence quark prediction of 3. 

The Ds/Fs ratio in VQCD is 

Ds/Fs{C.I.) = (95) 

From Rs = 3.086(19) for the smallest quark mass (k = 0.162), we obtain Ds/Fs = 
—0.021(4) which is close to zero as in the valence quark picture (Table 2) and differs 
from the lattice QCD calculation of - 0.58(18) (Sea + C. I.) and -0.31(11) (C. I.) (see 
Table 2) by a large margin. 
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6.3 Neutron to Proton Magnetic Moment Ratio 



We next come to the neutron to proton magnetic moment ratio The results 

of VQCD (indicated as A in Fig. 23) are very close to the well-known SU (6) value of 
- 2/3 (the result for the smallest quark mass case is - 0.662(22)), indirectly verifying 
the cloud effects of QCD (o for the C. I. in Fig. 23) which shows a 2.5 a departure 
from - 2/3 at the chiral limit. If there is any deviation of the VQCD from - 2/3, it 
should be due to the residual spin-spin interaction between the quarks in the baryon. 
Given the size of the error in our present results, we cannot make a definite conclusion 
on this aspect. 
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Figure 23: The magnetic moment ratio /i"/ /i^ between neutron and proton in VQCD 
and QCD are plotted against the dimensionless quark mass m^a from the strange to 
the charm region. A indicates the VQCD case and o/» indicates the C. I./ Sea -|- C. 
I. in QCD case. The dashed line is the SU{Q) prediction of - 2/3. 



7 Form Factors 

Now, we revisit the form factors of the nucleon. We calculated the isovector-axial form 
factor g\{q^)., the isoscalar-scalar form factor g%{q^)., the proton electric form factor 
G^E^'f) the isovector pseudoscalar form factor g%{q^) in VQCD at k = 0.162, 
which corresponds to the quark mass of ~ 120 MeV. They are plotted in Fig. 24 as 
a function of —q^. For comparison, we also plot in Fig. 25 the corresponding form 
factors from QCD at k = 0.154, which is about the same quark mass as in the VQCD 
case. 
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We see that although these form factors in VQCD are still different among them- 
selves, the differences are relatively smaller compared to those in QCD first of all, 
and secondly they are overall harder (except for gs{(l'^)) , i-e. they fall off slower than 
the corresponding ones in QCD. The most dramatic change is the pseudoscalar form 
factor where the size as determined by 

(r=) ^ -6^1.^0. (96) 

is reduced by about a factor of two. This is consistent with the pseudoscalar meson- 
dominance picture in Fig. 11, where the pseudoscalar form factor in QCD is domi- 
nated by the pion which in turn couples to the baryon through the nNN vertex. Yet, 
this meson 'cloud' is removed in VQCD by prohibiting the pair creation. In this case, 
the current couples directly to the quarks and consequently the (r^) of the hadron 
becomes smaller. To a lesser extent, similar situations happen in the vector and axial 
channels. This is again an indirect way of visualizing the effects of the meson clouds 
in the C. I. of QCD. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
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Figure 24: The isovector axial form factor g\{q^)i the isoscalar scalar form factor 
g%{q^)i the proton electric form factor G^^{q^) and the isovector pseudoscalar form 
factor gp{q^) in VQCD at k = 0.162 which corresponds to the quark mass of ~ 120 
MeV are plotted in terms of — g^. They are normalized at = to 1 in order to 
compare their dependence. 
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Figure 25: For comparison, the same 9\{q^), 9^s{'f)^ ^^(5^)5 ^^^d 9p{q'^) in QCD at 
K = 0.154 which is at about the same quark mass, i.e. ~ 120 MeV are plotted in 
terms of — o^. They are also normalized to 1 at (/^ = 0. 
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Figure 26: The neutron electric form factor G^{q^) for VQCD (•) at k = 0.162 is 
compared with the QCD result (o) at k = 0.154. 
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We also plot the neutron electric form factor G^{q^) for VQCD at k = 0.162 and 
its counter part (C. I. QCD at k = 0.154) in Fig. 26. We see that these two results 
are comparable in size and indicate that there are still some spin-spin correlation 
between the quarks in VQCD which breaks the SU{Q) symmetry. 



7.1 Goldberger-Treiman Relation and Vector Dominance 

There are several interesting aspects to observe in VQCD. Since the axial Ward iden- 
tities in Eqs. (^2|) and (|43|) are associated with the current involving both the u and 
V fields, they are applicable only to meson states with the creation and annihilation 
of quark- ant iquark pairs. Thus, the identities are useful in addressing the relation of 
the 'pion' mass and decay constant with the quark mass as PCAC is in QCD. On 
the other hand, it does not apply to baryons where only the quarks are involved. For 
example, the pseudoscalar current matrix element between the nucleon states does 
not have the pion pole as evidenced in Fig. 24. Consequently, there is no Goldberger- 
Treiman relation in VQCD. Conversely, the conserved vector current in Eq. ( p6D 
between the baryons and meson states leads to separately conserved quark and an- 
tiquark numbers. This entails the three-point function calculation as illustrated in 
Fig. 5(a). Yet, it does not apply to situations involving quark-antiquark creations 
or annihilations. For example, the vector meson decay constant (in the continuum) 
is determined by the matrix element of the vector current between the vector meson 
state and the vacuum 

mMk)) = ^ , (97) 

Jv 

where is the polarization vector, rriy is the vector-meson mass, and fy is the decay 
constant. Due to the vector current conservation 9^V^ = 0, the transversality condi- 
tion fc^ ■ = is derived in QCD. However in VQCD, the CVC relation in Eq. ( p6D 
does not apply to Eq. (0) because the conserved vector current in Eq. (^) does 
not have the pair annihilation term W7^u. Similarly, there is no vector dominance in 
the pion and nucleon EM form factors. As discussed in the preceeding section, there 
should be no meson dominance in form factors in VQCD. 

More generally, one can say that there is neither crossing relation, dispersion 
relation, nor unitarity in VQCD. But these features, or the lack of them, are shared 
by the valence quark model that we are trying to emulate. 



7.2 Quark Condensate and Symmetry Breaking 

The quark condensate in QCD is the order parameter of chiral symmetry break- 
ing. In VQCD, there are two quark condensates (uu) and (vv). Although VQCD 
eliminates the quark loops in the time direction, it still permits quark loops in the 
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spatial direction. One expects that 



\{uu)\ = \{vv)\ < (98) 

due to the fact that (uu) and (vv) have null measure in the time direction. Unfortu- 
nately there is an subtraction in the lattice calculation of quark condensate due to a 
contact term for the Wilson action |78[ . Short of a precise control of this subtraction, 



we are not able to have a fair comparison between {uu) and 

Nevertheless, (uu) and (vv) are nonzero which indicates that, to the extent that 
they serve as the order parameter of axial symmetry breaking as suggested by the 
existence of the Goldstone bosons and non-vanishing f^^ in the last section, the sym- 
metry breaking seems to be weaker than in QCD. It is shown in a Schwinger-Dyson 
equation study recently that the pseudoscalar meson grows either linearly or as 



the square root of the quark mass depending on whether it is large or small compared 
to a scale set by the quark condensate The linear dependence we see between 

the pion mass and the quark mass in Fig. 19 may mean that the quark masses we 
are calculating are still larger than the scale set by the quark condensate (uu) and 
(vv) and the quadratic pion mass dependence of the quark mass may yet to set in 
at smaller quark masses. Either way, the nonzero {uu) and {vv) supports the spon- 
taneous axial symmetry breaking scenario with two Goldstone pions as alluded to in 



section 5.4. 



8 Hadron Spectroscopy 

To further explore the consequences of the valence approximation, we study the 
hadron masses. Since hadron masses entail calculations of two-point functions, the 
sea quarks do not appear explicitly as they do in three-point functions (see Fig. 5(b)). 
The only exception is the flavor- singlet meson where the D. I. (Fig. 4(b)) is part of 
the meson propagator. The implicit sea quark effects in the loops which manifest 
themselves through the fermion determinant are known to affect the scaling [^, the 
topological susceptibility, phase transition, the r]' mass, and the slope of the hadron 
mass with respect to the quark mass |^T], 5U|. 



Here, we shall concentrate on the effects of the cloud quarks on hadron masses 
which are practically unknown. We first plot in Fig. 27 the masses of A,N,p, and 
TT as a function of the quark mass ma = \n{AKc/ k — 3) on our lattice with quenched 
approximation. We see that the hyper-fine splittings between the A and N, and the 
p and TT grow when the quark mass approaches the chiral limit as expected. 

In the infinite volume and continuum limits, it is found ||7^ that using nip to set 



the scale, the K, $ mesons and the octet and decuplet baryon masses are all within 
about 6% of the experimental results. 
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Figure 27: The dimensionless A,N,p, and vr masses in quenched QCD are plotted 
as a function of the quark mass m^a = ln(4Kc//t — 3). The pion mass is propor- 
tional to y/rriqa, while the others are extrapolated to the chiral limit with a linear m 
dependence. 

Next we plot in Fig. 28 the masses of A, A^, p, and vr as a function of the quark 
mass rriqa = ln{4K,c/i^ — 3) {nc = 0.1649 in this case) on the same quenched lattice. 
It is to our surprise that the truncation of the Z-graphs appears to have a dramatic 
effect on all these meson and baryon masses. 
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Figure 28: The dimensionless A, A^, p, and tt masses in VQCD are plotted as a function 
of the quark mass m^a = \n{4Kc/K — 3). All the masses are extrapolated to the zero 
quark mass limit with a linear niq dependence. 
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First of all, we notice that the A and the nucleon are on top of each other within 
errors all the way down to the smallest quark mass around the strange quark range. 
Thus, the hyper-fine splitting is largely gone in VQCD. This is true also between the 
p and vr. Extrapolating to the zero quark mass limit, the p mass m^a is 0.054(8). 
With fl-i = 1.75 GeV, nip = 95(14) MeV in VQCD. This is a factor of 6.4 smaller 
than that in the quenched approximation which gives rripa = 0.343(6). Secondly, we 
see that A, N, and p all dropped their masses greatly compared to those in QCD 
(Fig. 27). The extrapolated m/^a = 0.102(14), mAro = 0.074(11) at the zero quark 
mass limit. They are much smaller than their corresponding values of 0.638(41) and 
0.536(13) in the quenched QCD calculation in Fig 27. Furthermore, the hyper-fine 
splitting between A and N is now 49(7) MeV which is ~ 3.7 times smaller than our 
quenched result of 179(11) MeV Q. 

For a more direct comparison to see how the nucleon and A masses drop, we plot 
the nucleon and A masses in VQCD and quenched QCD on the same figure. We see 
in Fig. 29 that the nucleon moves down from VQCD to quenched QCD case by about 
a constant amount ~ 0.4 which is about 700 MeV. Similarly, the A mass drops by 
the same amount for heavier quarks. For quarks around the strange, it drops further 
to meet the nucleon. We also plot the p and vr masses in VQCD and quenched QCD 
on Fig. 30. Analogous to the N - A situation, the vector meson drops by about an 
equal amount ~ 0.31 or 537 MeV; whereas, the pseudoscalar meson drops about 0.22 
or 380 MeV in the strange region and approaches zero in both the quenched QCD 
and VQCD cases. 



''^Our quenched result is smaller than the experimental A-N splitting of 298 MeV mainly due to 
the fact that our results are not in the infinite volume and continuum limits. It is shown that when 
these limits arc taken, the octet and dccuplct baryons are within 6% of the experimental values p7|. 
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Figure 29: The dimensionless N and A masses from QCD are compared with those 
in VQCD function of the quark mass. 
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Figure 30: The dimensionless vr and p masses from QCD are compared with those in 
VQCD function of the quark mass. 
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Shown in Fig. 31 are the ai and ao mesons calculated in quenched QCD and 
VQCD. We see that both mesons come down in mass from QCD to VQCD by a large 
amount, ai appears to be degenerate with vr and p in VQCD over the range of the 
quark mass in Fig. 31 and Fig. 28. However, we cannot be certain about this point, 
especially in view of the fact that the errors on oi for the three lightest quarks are 
quite large, oq, on the other hand, seems to be higher than the pion in this range of 
the quark mass. 
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Figure 31: The dimensionless ai and masses from QCD are compared with those 
in VQCD function of the quark mass. 
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8.1 Tip correlation 



There is a concern about the near degeneracy of N and A and p and tt. Since VQCD 
is not Lorentz invariant, is it possible that the interpolation field M7iW intended for 
generating p may project onto a vr state as well and causes the vector correlation 
function {v'-fju{t)u'-fiv{0)) have the same asymptotic fall off as that of the pseudoscalar 
correlation function {v'jlu(t)u'j5v{0)) ? To check this possibility, we calculate the 
cross correlation functions among the vector (V), pseudoscalar (P), axial (A), and 
scalar (S) interpolation fields. We see in Fig. 32 the cross correlation functions 
{V(t)P{0)) , {A(t)P{0)), and (S'(t)P(O)) have essentially zero overlap over the range 
of the Euclidean time t. Also shown in Fig. 33 are the cross correlation functions 
{A{t)V{0)), {S{t)V{0)), and {P{t)V{0)). Again the overlaps are zero and the signals 
are cleaner. 
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Figure 32: The cross correlation functions between the pseudoscalar (P) and the 
vector (V), the axial (A) and the scalar (S) channels are shown as a function of the 
time separation ta. 



For comparison, we also plot the correlation function in the pseudoscalar and 
vector channels (P(t)P(O)) and {V{t)V{0)) for both the quenched and the VQCD 
cases in Fig. 34. We see that they are indeed much larger than the cross correlation 
functions in Figs. 32 and 33. 
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Figure 33: The cross correlation functions between the vector (V) and the axial (A), 
the scalar (S), and the pseudoscalar (P) channels are shown as a function of the time 
separation ta. 
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Figure 34: (a) The direct correlation functions (P(t)P(O)) for both the quenched and 
VQCD as a function of the time separation ta. (b) The same as in (a) for {V{t)V{0)). 
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8.2 Origin of Hyper-fine Splitting 



We see that the hyper-fine sphtting between the A and nucleon has largely disap- 
peared in the light quark sector when we remove the cloud quark and antiqurk in the 
form of the Z-graphs. This is rather mysterious in that according to the usual lore, the 
hyper-fine splitting is primarily due to the color-magnetic coupling induced spin-spin 



interaction between the quarks p, . This color-magnetic coupling is related to the 
spatial motion of the quarks which should not be affected by the truncation of the 
Z-graphs which only constraints the quark motion in the time direction. Indeed, this 
color-magnetic coupling is explicitly shown as <? ■ -B in the Pauli spinor representation 
of VQCD in Eqs. (P^ ) and (^Tj). Furthermore, it is this a ■ B term which is fully 
responsible for the hyperfine splitting between T and ?7fe in the heavy quark system. 
The latter is proven by the lattice QCD calculation with the non-relativistic QCD 
action containing such a term in the form of a ■ B/2Mi, 0, ^, |^. 

This raises a question as to how effective the color spin interaction is as far as 
the hyper-fine splitting is concerned in the light hadron spectroscopy. The same 
question has been raised by Glozman and Riska ETI, R^. Upon studying the negative 



parity and positive parity excitations of the N, A and A spectra, they found that the 
reverse ordering of the positive and negative parity resonances of N and particularly 
A from those of the A cannot be accommodated with the color-spin structure of the 
pairwise interaction ■ Xjffi ■ aj, rather it is consistent with the fiavor-spin structure 
Af ■ Ajaj ■ cTj. This is so because fiavor-spin structure of A is different from that 
of N and A. Interpreting this as due to Goldstone boson exchange, they can fit 
the low-lying baryon spectrum with a confinement potential in addition and also the 
magnetic moments of the baryon octets by taking into account of the meson exchange 
currents iS^ 



A similar problem was encountered in searching for scalar diquark clustering in 
lattice hadron form factors [Q]. Significant scalar diquark clustering is predicted in 
quark models which rely on the hyperfine interaction of the OGEP to split the N and 
A. While significant mass splitting is seen in the lattice simulations of Ref. P^ , PB | 
there is no evidence of scalar diquark clustering. This results leads one to look for 
other sources of hyperfine splitting that do not necessarily lead to clustering in the 
wave function, such as meson exchange [Q]. 

Furthermore, it is well known from the light baryon spectrum that the spin-orbit 
interaction is much weaker [^, ^ |^, |I0|, |^ than that which accompanies the spin- 



spin interaction in the one-gluon-exchange picture 0. This is problematic for the 
gluon-exchange picture if it is to explain both the heavy quarkonia which require 
the spin-orbit interaction and the light baryons which require a much weaker one. 
However, Goldstone boson exchange does not have the spin-orbit interaction between 
the light quarks and hence has no problem in this regard. 

This Goldstone boson exchange picture appears to be quite consistent with what 
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we find in VQCD. The flavor- nonsinglet meson exchange between the quarks is rep- 
resented by the Z-graph depicted in Fig. 32. Since all the Z-graphs are removed 
in VQCD, there will be no meson exchanges between the quarks as a result. This 
can explain why the hyper-fine splitting between A and nucleon is greatly reduced 
in VQCD (Fig. 29). But this does not answer the question as to why the color- 
magnetic coupling induced spin-spin interaction is not as effective in light baryons as 
in heavy quarkonium. While we don't have a strong evidence for it, we note that one 
aspect of the light quark may contribute to the difference. Unlike the heavy quarks, 
the propagator of the light quarks in the background gauge field can fluctuate into 
color-singlet meson clouds leading to meson dominance in various form factors (see 
Fig. 11). The range of fluctuation depends on the Compton wavelength of the meson. 
The longest range is in the pseudoscalar channel with the pion cloud as evidenced 
in the softness of the pseudoscalar form factor of the nucleon (Figs. 10 and 15). By 
the same token, Goldstone boson exchange between the quarks in Fig. 35 can have a 
range commensurate with its Compton wavelength. On the other hand, the range of 
one gluon exchange is limited since the gluon is confined. If the range of Goldstone 
boson exchange is longer than the gluon confinement scale, the hyper-fine interac- 
tion from Goldstone boson exchange is likely to to more effective than that from the 
color-magnetic coupling. In other words, the light quarks in the baryon have larger 
separations than those between the quarks and antiquarks in heavy quarkonia and 
this could be the reason for the diminished color-magnetic coupling in light baryon 
due to the limited range of the confined gluons. 



Figure 35: The meson exchange between quarks in the baryon is depicted as a Z- 
graph. The anti-quark produced in the Z-graph forms a meson with another quark 
in this case which is 'exchanged' between the two quarks. 




i 
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8.3 Origin of Dynamical Quark Mass 



Another significant feature of the VQCD spectroscopy in Figs. 28, 29, 30, and 31 
is that all the hadron masses drop substantially from their counterparts in QCD 
(including pion at finite quark mass). For example, the nucleon moves down from 
940 MeV in QCD (we used this to fix the scale) to 130(19) MeV in VQCD; A drops 
from 1117(72) MeV to 179(25) MeV; and p drops from 600(11) MeV to 95 (14) MeV. It 
is well known that chiral symmetry breaking leads to a dynamical quark mass related 
to the quark condensate in addition to the existence of Goldstone bosons. This 



can be seen from 

{^^) = {^l'^r + ^r^l) (99) 

which mixes the left- and right-handed quarks and has the effect of a dynamical mass 
clS 3i result of the chiral symmetry breaking. 

To the extent that we can interpret the VQCD result as due to the drop of dy- 
namical or constituent quark mass, we can draw the following conclusions: 

1. It is usually assumed in valence quark model that constituent quark mass arises 
from the dressing of the glue and the sea quark-antiquark pairs. Since the hard 
glue dressing in VQCD is expected to be the same as in QCD, it is not likely 
to be responsible for the dropping of hadron masses in VQCD. Furthermore, 



the quenched lattice calculations can reproduce the p, -ftT, $ mesons and 
the octet and decuplet baryon masses to within about 6% of the experimental 
results. This is an indication that the quark loops which generate sea quark- 
antiquark are not the primary source for hadron masses either. Here we see 
from our lattice calculation of VQCD that the dynamical quark mass actually 
arises from the 'dressing' of the cloud quarks — quark-antiquark pairs in the 
connected insertion. 

2. Since we conclude that the hyper-fine splitting in the light baryons is largely 
due to the Goldstone boson exchanges, one would expect the dynamical mass 
generation of the chiral symmetry breaking to be prominently realized in the 
baryon spectroscopy also. In the chiral symmetry models, the dynamical mass is 
generally generated through the a - the chiral partner of the pion. For example, 
in the linear sigma model, the dynamical mass is given [^, |91| as 



m,,„ = ^(^^)„^ (100) 

where mo- and g^qq and rriu are the sigma mass and its coupling to the quark. 
This is represented as the sigma-quark tadpole diagram as illustrated in Fig. 
33(a). Similar mechanism exists in the four-fermion Nambu-Jona-Lasinio model 
In QCD, the quark-line diagram which corresponds to the sigma-quark tadpole 
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in Fig. 36(a) would like something in Fig, 36(b) which inevitably involves cloud 
quarks and anti-quarks in the Z-graph. 



Interpreting the tadpole Z-graphs involving mesons in the scalar channel as the 
ones responsible for the dynamical quark mass generation, it is consistent with 
what we observe in VQCD. Dropping Z-graphs in VQCD which include these 
tadpoles diminishes the coupling to the quark condensate (^\E') and leads to 
the downturn of all the hadron masses from QCD. However, there is still a class 
of tadpole diagrams which survive. These are the spatial moving quark loops 
restricted within time slices. They may still couple to {uu) and {vv). But since 
these condensates in VQCD are much smaller than that in QCD (see Fig. 26), 
the dynamical mass is much smaller. This can explain why the masses of A, N, 
and p are small but non-zero in VQCD (Figs. 28, 29, and 30). 





a 



(b) 



Figure 36: (a) Sigma-quark tadpole diagram in the linear sigma model which is the 
mechanism for dynamical mass generation in this model, (b) The quark line diagram 
of the Z-graph which corresponds to the sigma-quark tadpole in (a). 



The interpretation we offer for the hyper-fine splitting and the dynamical quark 
mass is reminiscent of the chiral quark model of Manohar and Georgi [P^ on which 



the phenomenological studies of baryon masses |8T|, |82|, ^ and baryon structure 
are based. Arguing that the chiral symmetry breaking scale ^xSB is higher than the 
confinement scale Aqcd, they propose the relevant dynamical degrees of freedom to 
be the fundamental quarks, gluons and the Goldstone bosons in an effective theory 
at intermediate scales between A^^^ and Aqcd- What we observe in VQCD seems 
to suggest that the scale for the structure of baryons falls just in this range so that 
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the coupling to Goldstone bosons and dynamical mass generation are evident when 
QCD and VQCD are compared. 

There are other suggestions for the flavor-spin structure of the quark-quark inter- 
action. These are induced by instantons pSl |96| . It is known that the instantons give 



rise to chiral symmetry breaking and generate dynamical quark mass associated with 
{^"^f) 98]. The point-to-point hadronic correlation functions in the instanton 



liquid model have been verifled by lattice QCD calculation ||100|| and the role 
of instantons is revealed through cooling ||101|| . Although its direct connection to the 
cloud degree of freedom vis-a-vis VQCD is less transparent, the instanton picture, be- 
ing the root of chiral symmetry breaking, is expected to reproduce the consequences 
of the chiral quark model. 



9 Symmetry Breaking 

It is well known that the chiral symmetry SU(Nf)l x SU(Nf)r x Uv(1) of QCD is 
spontaneously broken to the diagonal SUv(Nf) x Uv(l). VQCD, as we have learned 
in this study, has a different symmetry breaking pattern. It starts out with the 
U(2Nf) symmetry (see Sec. |4.3|) with vector and axial symmetries in the particle- 
antiparticle space. Given our lattice simulation, which shows that the pseudoscalar 
mesons corresponding to the interpolation flelds u'j^v and v'j^u become massless at 
the zero quark mass limit, the pion decay constant is non-zero (it may actually 
diverge as l/rriT^), and the condensate {uu) and {vv) do not vanish. We take this as the 
evidence for spontaneous breaking of the axial symmetry in Eqs. ( ^2|) and (^). This 
then leads to a Uq(NF) x Uq(NF) symmetry, which is the vector symmetry for the 
quarks and antiquarks separately. Again by virtue of the lattice simulation, we find 
that SU{6) relation holds quite well for the g\/ Qa (or Fa/ Da), 9%/ 9% (or Ds/Fs), and 
IJ.n/lJ'p ratios. Furthermore, the nucleon and A, and p and vr are nearly degenerate. All 
these indicate that, although the SU{6) breaking color-magnetic coupling is present 



in the VQCD action (Eqs. (p9|), (0), and (|3lD), its effects are small. As a result, 
VQCD has an approximate higher symmetry, i.e. Uq(2NF) x Uq(2NF) where the 2 
represents the spin subgroup SU(2). This Uq(2NF) x Uq(2NF) is just the non-chiral 



U(6) X U(6) symmetry of Dashen and Gell-Mann [|102|| for Nf = 3 with quarks and 



antiquarks in the (6,1) and (1,6) representation respectively. It was proposed as 
a 'good' symmetry for "stationary (i.e. bound) and quasi-stationary (i.e. resonant) 
states of hadrons at rest" . It is interesting to note that after stripping off the sea and 
cloud quarks, we find VQCD possesses the same symmetry as in the valence quark 
model with spin degeneracy. 
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10 Conclusion, Analogy to Shell Model and Many 
Body Theory 



Instead of simulating QCD, we have mutilated it. The valence QCD theory we have 
constructed does not honor Lorentz invariance. It also violates unitarity, dispersion 
relation, and cross symmetry. But these are the attributes shared by the valence 
quark model which we set out to understand and our purpose of this study is to 
sort out the roles the various dynamical quark degrees of freedom play in different 
observables. This is much like the study of the brain |. One tries to correlate the 
dysfunction of certain part of the body with the damage of a specific part of the brain 
to infer its controlling mechanism. 

After defining the valence, the cloud, and the sea quarks from the hadronic tensor 
in deep inelastic scattering, we have been able to follow these degrees of freedom to 
three-point and two-point functions which are relevant to the quark model at low 
energies. Upon eliminating the cloud quarks in the connected insertion with the 
help of the VQCD action and the sea quarks in the disconnected insertion with the 
quenched VQCD calculation, we find from the ratios of g\l q'a, 9%/ g%i [Jin/ l^p and the 
masses of N, A,p, and tt that there is an approximate S't/(6) symmetry in VQCD 
which emerges from shaking off the "dressing" cloud and sea quark- ant iquark pairs. 
Its symmetry breaking pattern is distinct from that of QCD. We summarize the 
symmetry breaking pattern of QCD and VQCD in the following chart: 



We should point out that the Uq{2Np) x Uq{2Np) symmetry due to the spin degener- 
acy is only approximately true. We see that the ratio g\l g\ in Fig. 21 is not exactly 
3/5 and the neutron electric form factor Gg(g^) in Fig. 26 is not exactly zero. The 
degeneracies between A and N and between p and vr are not perfect either. These 
indicate that there are still some color-magnetic coupling induced spin effects. Nev- 
ertheless, it is a fairly good approximate symmetry. What we have demonstrated in 
this study is that QCD has this approximate symmetry in its valence approximation 
a la VQCD action. To our best knowledge, this is the connection between QCD and 
the valence quark model. 

The relation of the valence quark model and QCD actually is analogous to that 
between the shell model of nuclei and the many body theory. It is perhaps instructive 
to point out the parallel developments in the history of nuclear physics and hadron 
physics as far as the fermion dynamical degrees of freedom are concerned. We recall 
that the raison d'ete of the shell model consists of the pattern of energy levels, the 
spin and parity quantum numbers of nuclei, and the Schmidt lines for the magnetic 

■'■We thank T. Cohen for this analogy. 



QCD : 
VQCD : 



SU(Nf)l X SU(Nf)r X Uv(l) ^ SUv(Nf) x Uv(1) 

U(2Nf) ^ Uq(NF) X Uq(NF) Uq(2NF) X Uq(2NF) 
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moments of nuclei. Similar reasons, e.g. the mass pattern of baryons and mesons, 
SU(3) flavor symmetry, and the magnetic moments of proton and neutron lent their 
support for the existence of the quark model. Later experiments and theoretical 
developments in many body theory pointed out the inadequacies of the shell model 
and ideas like collectivity of the giant resonances ||103|| , pairing through the induced 



phonon-exchange interaction ||104|| , and core polarization for the magnetic moments 
or the Arima-Horie effect [|105|] are introduced. These involve the particle-hole degrees 
of freedom in the disconnected insertion which are the core polarization effects beyond 
the shell model. With the advent of QCD as the fundamental theory of quarks and 
gluons, similar ideas are introduced. For example, the resolution of the U(l) anomaly 
in terms of the topological susceptibility in the large Nc analysis by Witten and 
Veneziano [|19] is the schematic model ||103|| approach to generating the t]' mass by the 
collective coupling between quark loops. The concept of quark and gluon condensates 
are certainly related to pairing in the many body theory. Lack of appreciation for 
vacuum polarization due to the sea quarks for fiavor-singlet observables in the quark 
model has led to the "proton spin crisis" |]2D|, ^ and the nNa term puzzle |^ ^ pT| . 
The importance of Z-graphs for density dependence was pointed out for the effective 



nucleon-nucleon interaction ||106|| , and the higher density effects in the relativistic 
mean-field theory are shown largely to be due to the Z-graph with sigma meson 
exchanges ||107|| . The importance of the cloud quarks in hadrons through the Z- 



graphs are just beginning to be unraveled. The violation of the Gottfried sum rule 
leading to u{x) 7^ v{x) is shown to be due to the cloud antiquarks Furthermore, 
we have learned in the present study that the hyper-fine splitting in baryons and the 
dynamical quark mass are related to the cloud degree of freedom, which are probably 
the most surprising results of VQCD. 

The valence quark model, as we come to realize it today, is just like the shell 
model in nuclear physics. The U{6) x U{6) symmetry which comes with the valence 
quark model as the defining characteristic is not as good a symmetry as one tends to 
believe. Even with the supplement of SU{6) breaking one-gluon exchange, it does not 
capture the richness of the cloud degree of freedom in various form factors and matrix 
elements in the connected insertions. Moreover, the lack of sea degrees of freedom in 
the disconnected insertions is responsible for its overestimate of the fiavor-singlet (7^ 
by a factor of ~ 3 as well as its underestimate of the nNa term by a factor ~ 3 - 
4. What it lacks appears to be the spontaneous breaking chiral symmetry of QCD. 
This is exemplified in hadron spectroscopy where we find that the hyper-fine splitting 
between N and A and the dynamical quark mass are related to the cloud quark in 
the Z-graphs. 

One lesson we learned in this study is that the valence quark model is not nec- 
essarily a bad place to start building an effective theory of hadrons, provided one 
knows how to incorporate chiral symmetry and restores the cloud and sea degrees of 
freedom. Working in the intermediate scale between chiral symmetry breaking and 
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confinement which is appropriate for studying hadron structure and spectroscopy, one 
may start with the chiral quark model [p^. Integrating out the short-range part of the 



quark field is shown to lead to a very successful effective chiral theory of mesons [ 108 
One may extend this to the baryon sector with the quark coupling to the gluons and 
mesons . In this way, the cloud degrees of freedom will show up in the form factors 
and matrix elements via meson dominance and meson exchange currents. It can also 
give rise to the hyper-fine splitting and the dynamical quark mass. The meson loops 
on the quark lines, on the other hand, are responsible for the sea degrees of freedom. 
We will study this in the future. 
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